🚀 go-pugleaf

Mateuszzz pisze:
> Witam,
> 
> skoro z modułu GDAL używa większość programów GISowych, to może ktoś 
> zetknął się z następującym problemem.
> 
> Z nieznanej mi przyczyny nie mogę wykonać przeskalowania
> rastra w formacie ASCII (próbowałem przed konwersją zapisać go w postaci 
> GeoTiff - to również nie pomaga):
> 
> gdal_translate -outsize 5% 5% -of AAIGrid /home/costam.asc 
> /home/costam2.asc
> 
> Za każdym razem widzę błąd dotyczący nieprawidłowych wymiarów
> (0x0 TIFF). Co może być przyczyną?

A podejrzałeś tego grida w jakimś edytorze tekstowym? Wygląda na to, że 
  grid nie ma podanej wartości ncols i nrows (dwa pierwsze wpisy).

-- 
marcin

> A podejrzałeś tego grida w jakimś edytorze tekstowym? Wygląda na to, że
> grid nie ma podanej wartości ncols i nrows (dwa pierwsze wpisy).
>

GRID stworzyłem "samodzielnie" (napisałem skrypt, który wypisuje 
odpowiednie wartości do pliku tekstowego).
Jego (pliku .asc) zawartość to:

ncols 9
nrows 5
xllcorner 19.875
yllcorner 51.875
DX 0.25
DY 0.25
NODATA_value  -9999
1.9922307E-5, 1.894646E-5, 1.8956462E-5, 1.9361845E-5, 1.9511384E-5, 
1.9329693E-5, 1.9673384E-5, 1.8383538E-5, 1.8088462E-5
1.7483077E-5, 2.0421385E-5, 1.9258154E-5, 1.9281846E-5, 1.8062461E-5, 
1.7116461E-5, 1.5586462E-5, 1.8733692E-5, 1.9904E-5
1.781323E-5, 1.8908308E-5, 2.055877E-5, 1.998077E-5, 1.9613231E-5, 
1.8095077E-5, 1.7196462E-5, 1.5384308E-5, 1.757877E-5
1.9018922E-5, 1.895046E-5, 1.830677E-5, 2.0832154E-5, 2.051923E-5, 
2.0070462E-5, 1.8988E-5, 1.7819077E-5, 1.753277E-5
1.7792923E-5, 1.8765691E-5, 1.9184923E-5, 1.8667384E-5, 2.1397693E-5, 
2.1273538E-5, 2.0916308E-5, 1.9450616E-5, 1.8722461E-5

Mateuszzz pisze:
>> A podejrzałeś tego grida w jakimś edytorze tekstowym? Wygląda na to, że
>> grid nie ma podanej wartości ncols i nrows (dwa pierwsze wpisy).
>>
> 
> GRID stworzyłem "samodzielnie" (napisałem skrypt, który wypisuje 
> odpowiednie wartości do pliku tekstowego).
> Jego (pliku .asc) zawartość to:

Wygląda na ok, jedyne co może być problematyczne, to DX i DY (po 
skalowaniu przyjmują różne wartości), część programów tego nie 
obsługuje, bezpieczniej zastąpić jednym parametrem cellsize. A jakiej 
implementacji i wersji GDAL-a używasz? U mnie wtyczka do QGIS oparta na 
GDAL 1.7.3 wykonała bez problemu skalowanie Twojego pliku.

-- 
marcin

Mateuszzz pisze:

> GRID stworzyłem "samodzielnie" (napisałem skrypt, który wypisuje 
> odpowiednie wartości do pliku tekstowego).
> Jego (pliku .asc) zawartość to:
> 
> ncols 9
> nrows 5

Ok, już wiem, gdzie leży błąd. Jeżeli próbujesz przeskalować grida o 
wymiarach 9x5 o 5% to wynikowy będzie mieć rozmiary poniżej 1x1, a tak 
się nie da stąd zwraca Ci błąd. Albo większy grid albo mniejsze 
przeskalowanie.

-- 
marcin

> Ok, już wiem, gdzie leży błąd. Jeżeli próbujesz przeskalować grida o
> wymiarach 9x5 o 5% to wynikowy będzie mieć rozmiary poniżej 1x1, a tak
> się nie da stąd zwraca Ci błąd. Albo większy grid albo mniejsze
> przeskalowanie.
>

Hmmm...czyli najwidoczniej wpisuję złą wartość...
może opiszę całą sytuację....

chcę wykonać wykresy rozrzutu (oraz analizę statystyczną) danych 
znajdujących się w dwóch rastrach. Jeśli dobrze rozumuję, oba muszą mieć 
identyczne wymiary (myślę o współrzędnych rogów) i identyczną rozdzielczość.

W jaki sposób wykonać taką operację w QGISie?
Samo obcięcie nie jest kłopotliwe (Raster-->Kadruj, po czym można wpisać 
współrzędne rogów uzyskane przy pomocy narzędzia gdalinfo.

Nie wiem za bardzo jak "zmniejszyć" boki rastra (czy należy to 
zrealizować przez przeskalowanie?)
Kiedy wykonać taką operację - wydaje mi się, że logiczne byłoby 
zrobienie tego przed "kadrowaniem", jednak nawet w tej sytuacji może 
okazać się, że po obcięciu komórki jednego nie pokrywają się idealnie z 
bokami drugiego.

Wpadło mi do głowy użycie narzędzia gdal2xyz - zapisanie danych z rastra 
do pliku ze współrzędnymi oraz wartościami.
Następnym krokiem byłoby dodanie dodatkowych punktów (np. w 1/4 
odległości od środka rastra) z przepisaną wartością.

Tak dużą ilość punktów można byłoby "obciąć" do rozmiarów pierwszego z 
rastrów, po czym wykonać interpolację metodą najbliższego sąsiada.
Czy jest to dobry pomysł?
W jaki sposób szybko "zwielokrotnić" ilość punktów w rastrze?

> Ok, już wiem, gdzie leży błąd. Jeżeli próbujesz przeskalować grida o
> wymiarach 9x5 o 5% to wynikowy będzie mieć rozmiary poniżej 1x1, a tak
> się nie da stąd zwraca Ci błąd. Albo większy grid albo mniejsze
> przeskalowanie.
>
Dlaczego jest to niemożliwe?
Chciałbym uzyskać rastry o jeszcze mniejszym "oczku"
(teraz mam 0,25 piksela - chcę mieć mniejsze).

Czy możliwe jest wykonanie tej czynności?