如何单独lat经度点到3多边形?谁能帮我在这吗?

2视图(30天)

显示旧的评论

Bujee 2020年10月31日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/631989-how-to-separate-the-lat-lon-points-into-3-polygon-can-anyone-help-me-on-this-please

评论道: Bujee2020年10月31日

答:接受 KSSV

                         clc
                        
                         清晰的所有
                        
                         文件=负载(MATLAB \ 20190808 _4450_95104 \ ' D: \ \ str_CNVdataex.mat多边形的);
                        
                         fidout1 = fopen (MATLAB \ 20190808 _4450_95104 \ ' D: \ \ all.cnv多边形的,' w ');
                        
                         为k = 1:长度(file.str_CNV)
                        
                         on2 = file.str_CNV (k) .Ori.Year;
                        
                         sar2 = file.str_CNV (k) .Ori.Month;
                        
                         udur2 = file.str_CNV (k) .Ori.Day;
                        
                         tsag2 = file.str_CNV (k) .Ori.Hour;
                        
                         min2 = file.str_CNV (k) .Ori.Minute;
                        
                         sec2 = file.str_CNV (k) .Ori.Second;
                        
                         lat2 = str2double (file.str_CNV (k) .Ori.Lat);
                        
                         cns2 = file.str_CNV (k) .Ori.Cns;
                        
                         lon2 = str2double (file.str_CNV (k) .Ori.Lon);
                        
                         cew2 = file.str_CNV (k) .Ori.Cew;
                        
                         dep2 = str2double (file.str_CNV (k) .Ori.Dep);
                        
                         mag2 = file.str_CNV (k) .Ori.Mag;
                        
                         gap2 = str2double (file.str_CNV (k) .Ori.Gap);
                        
                         rms2 = file.str_CNV (k) .Ori.Rms;
                        
                         %如果lat2 < = 47.22 & & lat2 > = 46.70 & & lon2 < = 98.5 & & lon2 > = 97.00
                        
                         如果lat2 < = 49 & & lat2 > = 48 & & lon2 > = 95 & & lon2 < = 103 lat2 < = 47 & & lat2 > = 45.5 & & lon2 > = 95 & & lon2 < = 103
                        
                         %如果lat2 < = 48 & & lat2 > = 46 & & lon2 < = 102 & & lon2 > = 100如果lat2 = = 47.27 & & & & lat2 lon2 = = 98.24 = 98.00 = 47.00 & & lon2 = = & & & & lon2 lat2 = = 46.00 = = 101 & & & & lon2 lat2 = = 45.75 = = 100.76;
                        
                         %如果lat2 < = 47.27 & & lat2 > = 47.00 & & lon2 < = 98.24 & & lon2 > = 98.00 & & lat2 < = 46.00 & & lat2 > = 45.75 & & lon2 < = 101.00 & & lon2 > = 100.76
                        
                         流(fidout1' % s % s % s % s % s % s % 7.4 f % s % 8.4 f % s % 5.2 f % s % 3.0 f % s \ r \ n”on2, sar2 udur2、tsag2 min2, sec2, lat2, cns2, lon2, cew2, dep2, mag2, gap2, rms2);
                        
                         n = 0;
                        
                         为m = 1:长度(file.str_CNV (k) .ARR);
                        
                         sta = file.str_CNV (k) .ARR .STA (m);
                        
                         阶段= file.str_CNV (k) .ARR .Phase (m);
                        
                         魏= file.str_CNV (k) .ARR .Wgt (m);
                        
                         ttime = file.str_CNV (k) .ARR .TTime (m);
                        
                         流(fidout1“% s % s % s % s”sta阶段,魏,ttime);
                        
                         n = n + 1;
                        
                         如果n = = 6
                        
                         流(fidout1“\ r \ n”);
                        
                         n = 0;
                        
                         结束
                        
                         如果m = =长度(file.str_CNV (k) .ARR) & & n < 6 & & n ~ = 0
                        
                         流(fidout1“\ r \ n”);
                        
                         结束
                        
                         结束
                        
                         流(fidout1“\ r \ n”);
                        
                         %结束
                        
                         %结束
                        
                         结束
                        
                         结束
                        
                         流(fidout1“\ r \ n”);%输入

这是我的脚本读取我的盒子和排序的数据或矩形形状的区域。现在我想分离这些数据(dataex。三种,其中包括数据)三个区域。rectanugular区域和数据包括在该地区应分为3部分由两个obliguity或美德。

49.5 n / 95 e 49.5 n / 103 e 49 103 n / 95 e 47 n / e

49 n / 95 e 47 n / 103 e 48 n / 95 e 45.5 n / 103 e

48 n / 95 E 45.5 n / 103 E 45 103 n / 95 E 45 n / E

MATLAB文件=负载(' D: \ \多边形\ str_CNVdataex.mat”);这是所有的数据

之后,我想两个创建数据来自北方,中心,分别和南部的部分。

fidout1 = fopen (' D: \ MATLAB \多边形\ north.txt ', ' w ');

% fidout1 = fopen (' D: \ MATLAB \多边形\ center.txt ', ' w ');

% fidout1 = fopen (' D: \ MATLAB \多边形\ south.txt ', ' w ');

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

KSSV 2020年10月31日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/631989-how-to-separate-the-lat-lon-points-into-3-polygon-can-anyone-help-me-on-this-please answer_529794

你有限制....... (x0, y0)和(x1, y1)。让你所需网格分辨率和使用 interp2 。

                              x = linspace (x0, x1);
                             
                              y = linspace (y0, y1);
                             
                              (X, Y) = meshgrid (X, Y);
                             
                              Z = interp2(经度、纬度数据,X, Y);

3评论
显示2年长的评论隐藏2年长的评论

Bujee 2020年10月31日

因为我不理解为什么不是因为你的帮助。我非常感谢你的回答。我不是完整的新手,但我不能理解的事情相关的多边形。你能完成我的脚本如果可能吗?

登录置评。

类别

数学和优化映射的工具箱地图显示

雷达映射的工具箱地图显示

找到更多的在地图显示在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!