Mapping a raster with latitude and longitude coordinates

Question

Peter Dorofy 2015-8-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/232693-mapping-a-raster-with-latitude-and-longitude-coordinates

编辑： Peter Dorofy 2015-8-13

I have a satellite image file (2034x1354) where I pulled out reflectance values,latitude,longitude coordinates. I am having a difficult time trying to display/map this image. I have done research online and still cannot figure it out.

Here are the results:

Ref = 2030x1354 single
latitude = 2030x1354 single
longitude = 2030x1354 single

and reshaped to Vectors:

vRef = 1x2748620 single
vLat = 1x2748620 single 
vLon = 1x2748620 single

I have tried:

[latgrid,longrid] = meshgrid(latitude,longitude)
geoshow = (latgrid,longrid,Ref,'DisplayType','texturemap')

Not surprisingly, I get an error "Requested 2748620x2748620 (28144.2GB) array exceeds maximum..."

Any ideas?

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2015-8-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/232693-mapping-a-raster-with-latitude-and-longitude-coordinates#answer_188540

在 MATLAB Online 中打开

I would expect

geoshow(latitude, longitude, Ref, 'DisplayType', 'texturemap')

2 个评论
显示无隐藏无

Rob Comer 2015-8-9

I agree. There's no need to call meshgrid, because your 2D latitude-longitude arrays already constitute a mesh. Skipping this step will eliminate the latgrid and longrid arrays. That's good, because even if they fit in memory, these arrays would be both enormous and ill-formed.

Peter Dorofy 2015-8-13

编辑：Peter Dorofy 2015-8-13

Thanks everyone. Walter that worked. I actually did this before I posted, but it was not working. Probably because I did not save the script before executing.

请先登录，再进行评论。