如何從matlab中的shapefile只讀取多邊形的一個子集

我想從我的shapefile只讀取多邊形的子集。我正在使用'shaperead'來讀取shapefile，但我似乎無法只拉出我想要的多邊形。我知道我應該使用「選擇」對的值參數，但爲例網上沒有意義對我說：如何從matlab中的shapefile只讀取多邊形的一個子集

S = shaperead('concord_roads.shp','Selector',... 
     {@(v1,v2) (v1 >= 4) && (v2 >= 200),'CLASS','LENGTH'})

有一個在我的shapefile一個「站號」屬性，表明我的每個多邊形。我想要做的是能夠指定映射哪些多邊形（基於以前的聚類分析）。

來源

2016-07-04 new2matlab

wanted_stations = [...]; % as previously defined 

[S,C] = shaperead(filename,'selector',{@(v1) any(ismember(wanted_stations,v1)) ,'station_number'}); % assuming the attribute name is 'station_number'

基本上，對於每個多邊形，shaperead將屬性「station_number」的值賦給v1，然後評估匿名函數。 wanted_stations向量的值被合併到該函數中，如其定義的那樣。

來源

2016-07-04 21:07:19

非常感謝，這工作......雖然我必須在station_number周圍包括撇號和花括號以便它被識別。（即：wanted_stations = [{'25124'}]）...否則，形狀讀取功能只是繼續拉下所有東西。包含大括號使得數據類型成爲一個單元串，這導致了問題，因爲我想自動化這個過程，並且我的向量中的station_number類型是'cell'，但不是'cellstr'（我可以'似乎將它們轉換）。你有什麼建議嗎？ – new2matlab

看起來站號是作爲字符串存儲在shapefile中的，而不是整數，例如'25124'而不是25124.在這種情況下，如果您將該行修改爲any（ismember（wanted_stations，str2num（v1））），那麼您應該能夠將wanted_stations保留爲數字向量。 –

完美的工作！非常感謝。 – new2matlab

如何從matlab中的shapefile只讀取多邊形的一個子集

回答

相關問題