{"id":196,"date":"2021-12-10T08:47:06","date_gmt":"2021-12-10T00:47:06","guid":{"rendered":"http:\/\/purplepoolobservatory.com\/?p=196"},"modified":"2021-12-10T08:47:06","modified_gmt":"2021-12-10T00:47:06","slug":"%e5%a6%82%e4%bd%95%e7%94%a8matplotlib%e7%94%bbdem%e5%9b%be%e5%83%8f","status":"publish","type":"post","link":"https:\/\/purplepoolobservatory.com\/?p=196","title":{"rendered":"\u5982\u4f55\u7528matplotlib\u753bdem\u56fe\u50cf"},"content":{"rendered":"\n<pre class=\"wp-block-code\"><code>https:&#47;&#47;github.com\/rveciana\/geoexamples\nimport matplotlib.pyplot as plt\r\n\r\nfrom mpl_toolkits.basemap import Basemap\r\n\r\ndef hillshade(array, azimuth, angle_altitude):  # \u523b\u753b\u89c6\u89c93D\u6548\u679c\r\n\r\n    x, y = gradient(array)\r\n\r\n    slope = pi\/2. - arctan(sqrt(x*x + y*y))\r\n\r\n    aspect = arctan2(-x, y)\r\n\r\n    azimuthrad = azimuth*pi \/ 180.\r\n\r\n    altituderad = angle_altitude*pi \/ 180.\r\n\r\n    shaded = sin(altituderad) * sin(slope)\\\r\n\r\n     + cos(altituderad) * cos(slope)\\\r\n\r\n     * cos(azimuthrad - aspect)\r\n\r\n    return 255*(shaded + 1)\/2\r\n\r\nfig,ax=plt.subplots()\r\n\r\nm=Basemap(projection='merc',llcrnrlat=29,llcrnrlon=102,urcrnrlat=30,urcrnrlon=104,ax=ax)\r\n\r\n#ds = gdal.Open('\/home\/loong1\/data\/Asia_topo30.grd')\r\n\r\nds = gdal.Open('leshan.grd') #\u4e8b\u5148\u51c6\u5907\u597dgrd\u6570\u636e\u6587\u4ef6\r\n\r\nband = ds.GetRasterBand(1)\r\n\r\narr = band.ReadAsArray()#\u53ef\u4ee5\u4ecegrd\u6587\u4ef6\u91cc\u8bfb\u53d6\u90e8\u5206\u5185\u5bb9\uff0c\u4f46\u8fd9\u91cc\u8bfb\u53d6\u4e86\u5168\u90e8\r\n\r\nlon,lat=103,29.5\r\n\r\nx,y=m(lon,lat)\r\n\r\nhs_array = hillshade(arr,315, 45)\r\n\r\nm.scatter(x,y,marker='o',c='r')\r\n\r\n(x1,y1)=m(102,29)\r\n\r\n(x2,y2)=m(104,30)\r\n\r\nm.imshow(hs_array,cmap='Greys',origin='upper',extent=&#91;x1,x2,y1,y2]) # extent\u53c2\u6570\u6307\u5b9a\u56fe\u7247\u7684\u5750\u6807\u8303\u56f4\uff0c\u53e0\u52a0\u5230\u4e86basemap\u7ed9\u5b9a\u7684\u6295\u5f71\u5e95\u56fe\u4e0a\r\n\r\nm.drawparallels(arange(29,30,0.5),labels=&#91;0,1,1,0])\r\n\r\nm.drawmeridians(arange(102,104,1),labels=&#91;1,0,0,1])\r\n\r\nplt.show()<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"gallery","meta":{"footnotes":""},"categories":[2],"tags":[],"class_list":["post-196","post","type-post","status-publish","format-gallery","hentry","category-technology","post_format-post-format-gallery"],"_links":{"self":[{"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/posts\/196","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=196"}],"version-history":[{"count":1,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/posts\/196\/revisions"}],"predecessor-version":[{"id":197,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=\/wp\/v2\/posts\/196\/revisions\/197"}],"wp:attachment":[{"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=196"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=196"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/purplepoolobservatory.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=196"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}