使用rmaskout函数时报错

Hanchen

1. <div class="blockcode">#Maskout layer
   
2. lchina = shaperead('J - J - J')
   
3. cn_mask = lchina.shapes()
   
4. 
   
5. obsdir = r'D:'
   
6. obsfn = os.path.join(obsdir, 'SURF_CHN_MUL_HOR_N-2016120100.TXT')
   
7. 
   
8. obsf = DataFrame.read_table(obsfn)
   
9. lons = obsf['Lon'][:]
   
10. lats = obsf['Lat'][:]
    
11. odata = obsf['VIS'][:]
    
12. 
    
13. lons = DimArray(lons.values)
    
14. lats = DimArray(lats.values)
    
15. odata = DimArray(odata.values)
    
16. #之前使用micaps数据运行正常，查看数据类型为DimArray
    
17. #所以这里将数据转换为DimArray，但是还是会报错
    
18. #Maskout observation data by China border
    
19. obs, lon, lat = rmaskout(odata, lons, lats, cn_mask）

复制代码

报错如下:
>>> run script...
Traceback (most recent call last):
  File "C:\Users\18351\Desktop\MeteoInfo\code\test.py", line 18, in <module>
    obs, lon, lat = rmaskout(odata, lons, lats, cn_mask)
  File "C:\Users\18351\Desktop\MeteoInfo\MeteoInfo_3.3.3\MeteoInfo\pylib\mipylib\geolib\migeo.py", line 366, in rmaskout
    r = GeometryUtil.maskout_Remove(data.asarray(), x.asarray(), y.asarray(), mask)
        at org.meteoinfo.ndarray.ArrayString.getDouble(ArrayString.java:342)
        at org.meteoinfo.ndarray.IteratorFast.getDoubleNext(IteratorFast.java:80)
        at org.meteoinfo.geometry.geoprocess.GeometryUtil.maskout_Remove(GeometryUtil.java:555)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)
org.meteoinfo.ndarray.ForbiddenConversionException: org.meteoinfo.ndarray.ForbiddenConversionException


想请问各位大神这是什么原因呀？（使用的数据已经附上）

MeteoInfo

参考这里：http://www.meteothink.org/docs/m ... ame/read_table.html

Hanchen

MeteoInfo 发表于 2022-3-15 12:26
参考这里：http://www.meteothink.org/docs/meteoinfolab/dataframe/dataframe/read_table.html

已经解决了！谢谢老师！