- 积分
- 73
- 注册时间
- 2005-5-26
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-3-21 20:45:06
|
显示全部楼层
来自 湖北武汉
本帖最后由 alexqxp 于 2011-3-21 20:48 编辑
学习下pcolorm函数。
h = pcolorm(Z)
h = pcolorm(Z,npts)
h = pcolorm(lat,lon,Z)
h = pcolorm(lat,lon,Z,PropertyName,PropertyValue,...)
h = pcolorm(Z) projects the data grid Z on a graticule grid the size of Z between the latitude and longitude limits of the current map axes. The handle h of the displayed surface can be returned.
h = pcolorm(Z,npts) results in a graticule grid defined by npts, which is a two-element vector of the form [latitude-points longitude-points]. The default npts is [50 100].
h = pcolorm(lat,lon,Z) allows three other methods of defining the graticule grid. If lat and lon are matrices, they represent the actual graticule vertices as might be returned by meshgrat. If vectors, they are the representative coordinates of the rows and columns, respectively, of such a grid. If they are two-element vectors, they are treated as latitude and longitude limits, and a graticule mesh the size of the default npts is calculated.
h = pcolorm(lat,lon,Z,PropertyName,PropertyValue,...) allows the input of property name/property value pairs to control the surface object properties. Any property supported by the standard MATLAB function surface except XData, YData, and ZData can be altered in this manner.
RemarksThis function warps a data grid to a graticule mesh, which itself is projected according to the map axes property MapProjection. The fineness, or resolution, of this grid determines the quality of the projection and the speed of plotting it. There is no hard and fast rule for sufficient graticule resolution, but in general, cylindrical projections need fewer graticule points in the longitudinal direction than do complex curve-generating projections.
Examples
load topo
axesm miller
pcolorm(topo,[30 30])
demcmap(topo)
tightmap
根据例子来看,将海面高度数据替换topo,或许就可以将经纬度对应的海面高度映射到日本近海。随着时间高度发生变化,最后形成随时间变化的海啸图。当然,仅仅是分析猜测。欢迎楼主指正。 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|