MeteoInfoLab脚本：TXT站点数据绘制等值线填色图

MeteoInfo

1191246942 发表于 2020-7-11 17:00
弱弱的问下，想把边框上经纬度的字体变大点，应该用什么命令呢？

用 xaxis 和 yaxis 函数，里面有 tickfontsize 参数来调整坐标轴字体大小

段傻瓜瓜

111.21 30.74 35
111.31 30.74 47
111.25 30.64 51
111.40 30.67 63
111.34 30.61 70

diva211

按照楼主的脚本运行，提示以下信息，最后貌似提到了投影方式，不知道是啥原因的错误啊。。。是我的shp文件有问题吗
>>> run script...
D:/MI/data/o3.csv
Traceback (most recent call last):
  File "D:\MI\scripts\o3.py", line 16, in <module>
    bou1_layer = shaperead('d:/MeteoInfo_2.2.7/MeteoInfo/map/jinan.shp')
  File "D:\MeteoInfo_2.2.7\MeteoInfo\pylib\mipylib\geolib\migeo.py", line 56, in shaperead
    layer = MILayer(MapDataManage.readMapFile_ShapeFile(fn, encoding))
        at org.meteoinfo.projection.info.ProjectionInfo.esriStringToProj4Params(ProjectionInfo.java:408)

        at org.meteoinfo.projection.info.ProjectionInfo.factoryESRI(ProjectionInfo.java:138)

        at org.meteoinfo.data.mapdata.ShapeFileManage.loadProjFile(ShapeFileManage.java:718)

        at org.meteoinfo.data.mapdata.ShapeFileManage.loadShapeFile(ShapeFileManage.java:194)

        at org.meteoinfo.data.mapdata.MapDataManage.readMapFile_ShapeFile(MapDataManage.java:232)

        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.locationtech.proj4j.InvalidValueException: org.locationtech.proj4j.InvalidValueException: Unknown projection: Gauss_Kruger

diva211

diva211 发表于 2020-8-18 15:28
按照楼主的脚本运行，提示以下信息，最后貌似提到了投影方式，不知道是啥原因的错误啊。。。是我的shp文件 ...

已解决，是shp投影方式的问题

guiguixiong

                               
登录/注册后可看大图

来自苹果APP客户端

zqq

请问，我照着楼主的代码画图，为什么出来的图填色超出了地图边界呢？我的填色区域是整个画布

静言_GRMC

zqq 发表于 2021-3-19 21:29
请问，我照着楼主的代码画图，为什么出来的图填色超出了地图边界呢？我的填色区域是整个画布

gtemp,gx,gy = griddata((lon, lat), pm, xi=(x, y), method='idw', radius=0.8)
radius=0.8 调整一下这个

zqq

静言_GRMC 发表于 2021-10-9 10:56
gtemp,gx,gy = griddata((lon, lat), pm, xi=(x, y), method='idw', radius=0.8)
radius=0.8 调整一下这 ...

好的！谢谢老师！

ybsky

fn = 'C:/Users/ybsky/Desktop/porject1/rain.txt'
ncol = numasciicol(fn)
nrow = numasciirow(fn)
a = asciiread(fn,shape=(nrow,ncol))
lon = a[:,0]
lat = a[:,1]
t = a[:,2]
#To grid data
x = arange(104, 106.7, 0.2)
y = arange(32.5, 34.7, 0.2)
gtemp,gx,gy = griddata((lon, lat), t, xi=(x, y), method='idw', radius=0.1)
#Plot
axesm()
geoshow('ln',edgecolor='lightgray')

levs = [0, 5, 10, 15, 20, 25, 35, 40, 45,50]
cols = [(255,255,255),(0,255,0),(127,255,0),(255,255,0),(255,215,0),(255,128,0),(255,97,0), \
    (255,0,0),(176,23,31),(135,38,87),(255,0,255)]
#layer = contourfm(x, y, gtemp,20)  
layer = contourfm(x, y, gtemp,levs,colors=cols)
#slayer = scatterm(lon, lat,t,colors=['k'], size=10)
slayer = scatterm(lon, lat, t,levs,colors=cols, size=8)
masklayer(mlayer, [layer])
xlim(104, 106.7)
ylim(32.5, 34.7)
title(u'陇南暴雨频次',fontname=u'黑体',fontsize=20,bold=False,color='blue')   
#text(113.3, 23.9, u'2017年5月6日13时', fontname=u'黑体', fontsize=16)
colorbar(layer)

代码正常运行
34.26 104.18   7
33.81 105.12   8
32.90 104.81   6
33.72 105.54   17
33.76 105.80   29
33.86 105.69   22
33.99 105.34   12
33.98 105.31   5
32.78 105.38   48
32.84 104.99   6
34.29 104.15   7
33.34 105.67   34
32.95 104.67   12
32.98 105.25   25
33.30 105.46   28
33.43 104.82   10
32.98 105.08   9
34.03 105.30   6
这是数据，为啥出图了只有地图，没有数据

MeteoInfo

ybsky 发表于 2022-5-6 16:25
fn = 'C:/Users/ybsky/Desktop/porject1/rain.txt'
ncol = numasciicol(fn)
nrow = numasciirow(fn)

lon, lat 弄反了