Temat: Reklasyfikacja kilku rastrów wg jednego schematu

Za pomocą skryptu można wyciągnąć jedno pasmo (na raz dla wszystkich map), dla ArcGIS 10 będzie to wyglądać tak:

import arcpy, os
from arcpy import env

# parametry wejściowe, tu trzeba wpisac folder, w ktorym sa rastry oraz folder, do ktorego zapiszemy wynik

arcpy.env.workspace = "C:/..."
outworkspace = "C:/..."

# numer pasma rgb - edit: trzeba było podać w nawiasach
bandno = "1"

rasterList = arcpy.ListRasters("*", "")

for raster in rasterList:
inraster = raster + "/Band_" + bandno
outraster = outworkspace + "/" + raster + "_b" + bandno + ".tif"arcpy.CopyRaster_management(inraster,outraster,"DEFAULTS","0","9","","","8_BIT_UNSIGNED")

Jeżeli już będzie się mieć rastry zapisane w formacie 8 bitowym, to można je również za pomocą skryptu zreklasyfikować - ale sadzę, że to nie wyjdzie dobrze. Odcienie szarości będą zapewne różne dla kolejnych map.

rasterList = arcpy.ListRasters("*", "")

for raster in rasterList:

reclassField = "Value"
#przykladowa reklasyfikacja 0 na 1 i 255 na 0
remap = RemapValue([[0, 1], [255, 0],])

arcpy.CheckOutExtension("Spatial")

outReclassify = Reclassify(raster, reclassField, remap, "NODATA")

outraster = outworkspace + "/" + raster + "_reklas" + ".tif"

outReclassify.save(outraster)

Jeżeli jest jakiś prostszy sposób, to chętnie go poznam.

[edit] był mały błąd - numer pasma trzeba podać w nawiasie, teraz działa (u mnie)Marcin Szwagrzyk edytował(a) ten post dnia 07.02.12 o godzinie 20:11

Temat: Reklasyfikacja kilku rastrów wg jednego schematu

Jak wyglądają takie diapozytywy? Są rozdzielone na warstwy? Na czym polega przewaga takiego diapozytywu (jeśli chodzi o automatyczną wektoryzaję) w stosunku do zwykłego rastra mapy topograficznej?

Kiedyś tam wyciągliśmy z kolegą zabudowę z mapy topo (10k) i -pomijając dziury wycięte przez napisy i inne obiekty nachodzące na budynki- wyszła trochę koślawa, ale całkiem przyzwoita warstwa wektorowa.

Pomysł z tymi diapozytywami ciekawy, chętnie dowiem się więcej.

[EDIT]
Na rastrze w skali szarości (wyciągnięte jedno pasmo) budynki mają taką samą wartość piksela, jak poziomice i coś tam jeszcze, więc się ich w ten sposób nie wyodrębni z reszty. Trzeba by pewnie stworzyć raster 8 bitowy, z wykorzystaniem wszystkich pasm, i z tego wyselekcjonować zabudowę...Marcin Szwagrzyk edytował(a) ten post dnia 07.02.12 o godzinie 20:37

Temat: Reklasyfikacja kilku rastrów wg jednego schematu

Zależy, do czego ta warstwa jest potrzebna. Jak np do jakiejś wizualizacji w skali 10k, to wyciągnięcie kolorów z rastra, zamiana na wektor (wywalenie śmieci z automatu - np obiektów za małych na budynki czy też obiektów o absurdalnych wskaźnikach obwód/powierzchnia), w zupełności wystarczy.

Nadzorowana klasyfikacja i tak nie da jakiś super rezultatów - problemem będą pewnie napisy zachodzące na budynki (w miastach jest ich sporo) i inne obiekty przecinające budynki - np linie energetyczne. Chyba, że potraficie sobie z tym poradzić.

Temat: Reklasyfikacja kilku rastrów wg jednego schematu

W każdym razie, żeby pokazać, że się da (czy jest to opłacalne i sensowne, to inna sprawa).

Pierwszy obrazek - po wyciągnięciu pikseli o wartości odpowiadającej zabudowie mieszkaniowej (z mapy w skali szarości).

Drugi obrazek - po obróbce - usunięciu obiektów mniejszych niż 20 m2, zbuforowaniu o 2 m (bo obcięło trochę z brzegów) i zastosowaniu narzędzia simplify buildings

http://imageshack.us/f/195/zabudowag.jpg/Marcin Szwagrzyk edytował(a) ten post dnia 08.02.12 o godzinie 14:41