如何将lat, long转换为xy在地图上绘图?

3次浏览(过去30天)

显示旧的注释

Sabereh 2022年10月8日

0
链接

直接链接到这个问题

https://es.mathworks.com/matlabcentral/answers/1820055-how-to-convert-lat-long-to-x-y-for-plotting-on-the-map

评论道: 本杰明•汤普森 2022年10月10日

我有一个14080*13624 .jpg图像，这是国家地图。

我还知道地图上四个角的纬度和长度。

我有其他点的纬度和长度我想把它们分散在地图上。

我使用下面的代码，但结果不匹配点的真实位置。

                          Img = imread(“map.jpg”）;
                         
                          imshow (img);
                         
                          持有在
                         
                          %%设计点
                         
                          拉特= 32.778473;
                         
                          Lon = 51.648468;
                         
                          H = 400;
                         
                          % %中心
                         
                          Lat0 = 31.8978;
                         
                          Lon0 = 53.6133;
                         
                          H0 = 400;
                         
                          % %
                         
                          wgs84 = wgs84椭球;
                         
                          [xEast,yNorth,zUp] = geodetic2enu(lat,lon,h,lat0,lon0,h0,wgs84);
                         
                          Pix_size_x = 14080/2;
                         
                          Pix_size_y = 13824/2;
                         
                          Stepx = (8.2281e+05)/pix_size_x;
                         
                          step py = (9.2292e+05)/pix_size_y;
                         
                          x = pix_size_x + xEast/stepx
                         
                          y = pix_size_y - yNorth/stepy;
                         
                          散射(x, y)
                         
                          轴平等的

你能帮我解决这个问题吗?

5个评论
显示隐藏 4旧的评论

沃尔特·罗伯森 2022年10月9日

稍微清理了一下代码，但没有解决问题。

对我来说，为什么把像素除以2不是很明显。与参考点位于地图中心有关吗?

我建议您使用地图边界的坐标作为输入，并让它计算中心和跨度。

                                %%映射大小
                               
                                Long_covered = 8.2281e+05;
                               
                                Lat_covered = 9.2292e+05;
                               
                                %%设计点
                               
                                拉特= 32.778473;
                               
                                Lon = 51.648468;
                               
                                H = 400;
                               
                                % %中心
                               
                                Lat0 = 31.8978;
                               
                                Lon0 = 53.6133;
                               
                                H0 = 400;
                               
                                % %
                               
                                Img = imread(“map.jpg”）;
                               
                                imshow (img);
                               
                                持有在
                               
                                wgs84 = wgs84椭球;
                               
                                [xEast,yNorth,zUp] = geodetic2enu(lat,lon,h,lat0,lon0,h0,wgs84);
                               
                                Pix_size_x = size(img,2)/2;
                               
                                Pix_size_y = size(img,1)/2;
                               
                                Stepx = long_covered/pix_size_x;
                               
                                Stepy = lat_covered/pix_size_y;
                               
                                x = pix_size_x + xEast/stepx
                               
                                y = pix_size_y - yNorth/stepy;
                               
                                散射(x, y)
                               
                                轴平等的