程序=算法+數(shù)據(jù)結(jié)構(gòu)

計(jì)算機(jī)是一個(gè)固定的一個(gè)程序段，或稱其為一個(gè)子程序，它在可以實(shí)現(xiàn)固定運(yùn)算功能的同時(shí)，還帶有一個(gè)入口和一個(gè)出口，所謂的入口，就是函數(shù)所帶的各個(gè)參數(shù)，我們可以通過(guò)這個(gè)入口，把函數(shù)的參數(shù)值代入子程序，供計(jì)算機(jī)處理；所謂出口，就是指函數(shù)的函數(shù)值，在計(jì)算機(jī)求得之后，由此口帶回給調(diào)用它的程序。

效果那是極好的。

• 函數(shù)名：draw_shp

• 輸入?yún)?shù)：shapefile的路徑

• 輸出：直接顯示地圖。

• 需要用的到的包：

• shapefile：讀取shp文件

• plotly：繪圖

• geo_types：序列化幾何對(duì)象 編碼實(shí)現(xiàn)設(shè)計(jì)：###在Cargo.toml文件里面，導(dǎo)入需要的包：

//Cargo.toml
//后面的features特性，暫時(shí)不用去管，這是一種Rust特有的編譯特性
plotly = { version = "0.8.4", features = ["kaleido"] }
shapefile = {version = "0.5.0", features = ["geo-types"]}
geo-types = "0.7.12"

讀取一個(gè)shapefile，并且把幾何信息給獲取出來(lái)。

在Rust中，可以通過(guò)shapefile包來(lái)讀取shapefile，實(shí)現(xiàn)如下：

let shp = shapefile::read_as::<_,
        shapefile::Polygon, shapefile::dbase::Record>(
            "./data/shp/北京行政區(qū)劃.shp",
    ).expect(&format!("Could not open polygon-shapefile: './data/shp/北京行政區(qū)劃.shp'"));

接下去，需要把里面的geometry信息給取出來(lái)：

//定義一個(gè)集合，通過(guò)文件迭代器，把geometry部分轉(zhuǎn)換成polygon，然后加入到這個(gè)幾何里面去。
let mut polygons:Vec<Polygon> = Vec::new();
for (polygon, polygon_record) in shp {
    let geo_mpolygon: geo_types::MultiPolygon<f64> = polygon.into();
    for poly in geo_mpolygon.iter(){
        polygons.push(poly.to_owned());
    }
}

把這個(gè)polygon集合，繪制到plotly上去

首先plotly繪制幾何圖形，是按照坐標(biāo)來(lái)的，一系列坐標(biāo)組成一個(gè)繪圖元素，代碼如下：

• 注1：plotly的地圖繪制用的是mapbox的api，所以是先維度lat，再經(jīng)度lon，得反過(guò)來(lái)。

• 注2：trace是plotly繪圖的基本元素，這里每個(gè)幾何要素（如一個(gè)面），就可以構(gòu)建一個(gè)trace，如果可以設(shè)置為一個(gè)顏色，也可以繪制為不同的顏色。

let mut trace_vec = Vec::new();
    for ps in polygons{
        let mut lon:Vec<f64> = Vec::new();
        let mut lat:Vec<f64> = Vec::new();
        for p in ps.exterior(){
            lon.push(p.x);
            lat.push(p.y);
        }
        let trace = ScatterMapbox::new(lat, lon).mode(Mode::None)
        .fill(plotly::scatter_mapbox::Fill::ToSelf)
        .fill_color(Rgba::new(0,0,255,0.5));
        trace_vec.push(trace);
    }

獲得基本地圖的配置

• 可以看見(jiàn)這里用的是mapbox的地圖，默認(rèn)的風(fēng)格為無(wú)背景的白板，默認(rèn)的中心位置是東經(jīng)116.3,北緯39.9，默認(rèn)地圖放大等級(jí)是9級(jí)。這些信息都是地圖的初始化的默認(rèn)配置。

let layout = Layout::new()
    .drag_mode(DragMode::Zoom)
    .margin(Margin::new().top(10).left(10).bottom(10).right(10))
    .width(1024)
    .height(700)
    .mapbox(
        Mapbox::new()
            .style(MapboxStyle::WhiteBg)
            .center(Center::new(39.9, 116.3))
            .zoom(9),
    )

最后繪制顯示地圖

let mut plot = Plot::new();
plot.set_layout(layout);
for t in trace_vec.iter(){
    plot.add_trace(t.to_owned());
}
plot.show();

全部代碼，放到一個(gè)function里面，如下所示：

然后寫個(gè)測(cè)試方法運(yùn)行一下：

運(yùn)行結(jié)果：