r.sun QGIS (GRASS)

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 11:09

list писал(а): Пробовал r.out.gdal, r.out.ascii, r.out.xyz, r.out.arc. Пробовал в QGIS перевести на другой формат (Gtiff, hdr.adf...), все равно потом АРК ГИС его не видит.
И что, ArcGIS и .tif (геотифф) не видит? Кстати, что значит "не видит"? Не может распознать СК/проекцию растра? Или не открывает растр (в чём я лично сильно сомневаюсь)?
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 11:18

Александр Мурый писал(а):
list писал(а): что значит "не видит"? Не может распознать СК/проекцию растра? Или не открывает растр (в чём я лично сильно сомневаюсь)?
Пусто там. Папка открывается а внутри пусто.

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 11:26

Если вы делали .tif через <r.out.gdal>, то приведите здесь вывод команды "gdalinfo имя_файла.tif".
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 12:09

Александр Мурый писал(а):Если вы делали .tif через <r.out.gdal>, то приведите здесь вывод команды "gdalinfo имя_файла.tif".
Вроде нашел корень проблем. Непосмотрел в вывод команд.
(Fri Nov 6 09:36:36 2015)
r.out.gdal input=dif_jar@PERMANENT output=/home/lis/Dokumenty/dp/smaz/a
export do datového typu GDAL: Float32
Checking GDAL data type and nodata value
Exporting to GDAL raster
ERROR 6: SetColorTable() only supported for Byte or UInt16
bands in TIFF format.
r.out.gdal hotovo.
(Fri Nov 6 09:36:39 2015) Příkaz ukončen (2 sec)
В чем разница между Float32, Byte или UInt16 если в двух словах сказать.
Float значит что можно после запятой писать числа (с плавающей десятинной). Что означает 32, возможно количество цифр после запятой.
Byte значит по-моему что только целые числа
UInt16 не знаю
Но что мне нужно, мне нужны десятинные значения. Тогда наверное UInt16

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 12:25

list писал(а): Вроде нашел корень проблем. Непосмотрел в вывод команд.
r.out.gdal input=dif_jar@PERMANENT output=/home/lis/Dokumenty/dp/smaz/a
export do datového typu GDAL: Float32
Попробуйте флаг "-f" (Force raster export despite any warnings of data loss).
list писал(а): В чем разница между Float32, Byte или UInt16 если в двух словах сказать.
Float значит что можно после запятой писать числа (с плавающей десятинной). Что означает 32, возможно количество цифр после запятой.
Byte значит по-моему что только целые числа
UInt16 не знаю
Не думаю, что у вас отключили гугл. UInt16 — 16-разрядное целое число без знака.
list писал(а): Но что мне нужно, мне нужны десятинные значения. Тогда наверное UInt16
Нужен Float32.
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 12:55

Александр Мурый писал(а):
list писал(а): Попробуйте флаг "-f" (Force raster export despite any warnings of data loss).
От ERROR избавился (-c), но в АРК ГИСе все по-прежнему. В папке пусто. А что интересно обратно в ГРАСС (в локацию) экспортированные растры не начитываются.

(Fri Nov 6 10:28:59 2015)
r.out.gdal -f input=dif_jar@PERMANENT type=Float32 output=/home/lis/Dokumenty/dp/smaz/d
export do datového typu GDAL: Float32
Checking GDAL data type and nodata value
Exporting to GDAL raster
ERROR 6: SetColorTable() only supported for Byte or UInt16
bands in TIFF format.
r.out.gdal hotovo.
(Fri Nov 6 10:29:02 2015) Příkaz ukončen (3 sec)
(Fri Nov 6 10:29:37 2015)
r.out.gdal -c -f input=dif_jar@PERMANENT type=Float32 output=/home/lis/Dokumenty/dp/smaz/e
export do datového typu GDAL: Float32
Checking GDAL data type and nodata value
Exporting to GDAL raster
r.out.gdal hotovo.
(Fri Nov 6 10:29:40 2015) Příkaz ukončen (2 sec)
(Fri Nov 6 10:36:46 2015)
r.out.gdal -c input=dif_jar@PERMANENT type=Float32 output=/home/lis/Dokumenty/dp/smaz/f
export do datového typu GDAL: Float32
Checking GDAL data type and nodata value
Exporting to GDAL raster
r.out.gdal hotovo.
(Fri Nov 6 10:36:49 2015) Příkaz ukončen (2 sec)

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 13:01

Приведите вывод команды "r.info dif_jar".
Вы задаёте вычислительный регион по растру перед экспортом?

Код: Выделить всё

g.region rast=dif_jar
Я правильно понял, что файла на выходе нет?
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 13:08

Александр Мурый писал(а): Я правильно понял, что файла на выходе нет?
По-моему да. В папке он есть но ни ГРАСС ни АРК ГИС его не видят

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 13:12

list писал(а): По-моему да. В папке он есть но ни ГРАСС ни АРК ГИС его не видят
Что значит "по-моему"? Кажущийся файл? :) Он либо есть, либо его нет. Посмотрите внимательно в каталог "/home/lis/Dokumenty/dp/smaz/e".
GRASS и не может его видеть (т.к. использует собственный формат, а не .tif).
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 13:19

Александр Мурый писал(а):
list писал(а): Посмотрите внимательно в каталог "/home/lis/Dokumenty/dp/smaz/e".
Там есть

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 13:22

Какой размер файла? Файл открывается, например, в QGIS?
Выполните в командной строке "gdalinfo /home/lis/Dokumenty/dp/smaz/e/ИМЯ_ВАШЕГО_ФАЙЛА.tif" и напишите сюда вывод команды.
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 13:25

Александр Мурый писал(а):Приведите вывод команды "r.info dif_jar".

Код: Выделить всё

 Layer:    dif_jar                        Date: Fri Oct 30 07:38:57 2015    |
 | Mapset:   PERMANENT                      Login of Creator: lis             |
 | Location: r.sun_test                                                       |
 | DataBase: /home/lis/Dokumenty/dp/pokus/location                            |
 | Title:     ( pri_jar )                                                     |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 255             |
 |   Data Type:    FCELL                                                      |
 |   Rows:         1500                                                       |
 |   Columns:      1875                                                       |
 |   Total Cells:  2812500                                                    |
 |        Projection: Krovak                                                  |
 |            N: -987250.03    S: -988750.03   Res:     1                     |
 |            E: -640625.03    W: -642500.03   Res:     1                     |
 |   Range of data:    min = 408.0876  max = 408.1714                         |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.sun                                                      |
 |                                                                            |
 |   Comments:                                                                |
 |     ----------------------------------------------------------------       |
 |     Day [1-365]:                              79                           |
 |     Solar constant (W/m^2):                   1367                         |
 |     Extraterrestrial irradiance (W/m^2):      1378.781555                  |
 |     Declination (rad):                        -0.002668                    |
 |     Latitude min-max(deg):                    50.6873 - 50.7027            |
 |     Sunrise time min-max (hr.):               6.01 - 6.01                  |
 |     Sunset time min-max (hr.):                17.99 - 17.99                |
 |     Time step (hr.):                          0.5000                       |
 |     Linke turbidity factor min-max:           4.6-4.6                      |
 |     Ground albedo:                            0.000                        |
 |     -----------------------------------------------------------------      |
 |                                                                            |
 |    r.sun -s elevin="dtm@PERMANENT" aspect=270 slope=0.0 linkein="jaro@P\   |
 |    ERMANENT" lin=4.6 alb=0.0 coefbh="beam@PERMANENT" coefdh="dif@PERMAN\   |
 |    ENT" beam_rad="pri_jar" insol_time="hod_jar" diff_rad="dif_jar" glob\   |
 |    _rad="glo_jar" day=79 step=0.5 dist=1.0 numpartitions=1   

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 13:49

Александр Мурый писал(а):Какой размер файла? Файл открывается, например, в QGIS?
QGIS тоже его не открывает. Размер файла 11,3 мегабайт.
Выполните в командной строке "gdalinfo /home/lis/Dokumenty/dp/smaz/e/ИМЯ_ВАШЕГО_ФАЙЛА.tif" и напишите сюда вывод команды.

Код: Выделить всё

 GRASS 6.4.3 (r.sun_test):~ > gdalinfo /home/lis/Dokumenty/dp/smaz/e
Driver: GTiff/GeoTIFF
Files: /home/lis/Dokumenty/dp/smaz/e
Size is 1875, 1500
Coordinate System is:
PROJCS["Krovak",
    GEOGCS["GCS_bessel",
        DATUM["D_Militar_Geographische_Institut",
            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
        PRIMEM["Ferro",-17.666666666668],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Krovak"],
    PARAMETER["latitude_of_center",49.5],
    PARAMETER["longitude_of_center",42.5],
    PARAMETER["azimuth",0],
    PARAMETER["pseudo_standard_parallel_1",0],
    PARAMETER["scale_factor",0.9999],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
Origin = (-642500.030000000027940,-987250.030000000027940)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( -642500.030, -987250.030) ( 33d23'17.71"E, 50d42' 2.44"N)
Lower Left  ( -642500.030, -988750.030) ( 33d23'26.81"E, 50d41'14.24"N)
Upper Right ( -640625.030, -987250.030) ( 33d24'52.59"E, 50d42' 9.65"N)
Lower Right ( -640625.030, -988750.030) ( 33d25' 1.66"E, 50d41'21.45"N)
Center      ( -641562.530, -988000.030) ( 33d24' 9.69"E, 50d41'41.95"N)
Band 1 Block=1875x1 Type=Float32, ColorInterp=Gray
  Metadata:
    COLOR_TABLE_RULE_RGB_0=4.080876e+02 4.081043e+02 255 255 0 0 255 0
    COLOR_TABLE_RULE_RGB_1=4.081043e+02 4.081211e+02 0 255 0 0 255 255
    COLOR_TABLE_RULE_RGB_2=4.081211e+02 4.081379e+02 0 255 255 0 0 255
    COLOR_TABLE_RULE_RGB_3=4.081379e+02 4.081547e+02 0 0 255 255 0 255
    COLOR_TABLE_RULE_RGB_4=4.081547e+02 4.081714e+02 255 0 255 255 0 0
    COLOR_TABLE_RULES_COUNT=5

Александр Мурый
Гуру
Сообщения: 5149
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: r.sun QGIS (GRASS)

Сообщение Александр Мурый » 06 ноя 2015, 13:52

Экспортируйте файл в файл, а не в каталог:

Код: Выделить всё

r.out.gdal -f input=dif_jar@PERMANENT type=Float32 output=/home/lis/Dokumenty/dp/smaz/e/ИМЯ_ФАЙЛА.tif
И попробуйте открыть где-нибудь ещё раз.
Редактор материалов, модератор форума

list
Интересующийся
Сообщения: 28
Зарегистрирован: 10 окт 2015, 01:53
Репутация: 0

Re: r.sun QGIS (GRASS)

Сообщение list » 06 ноя 2015, 14:15

Александр Мурый писал(а):Экспортируйте файл в файл, а не в каталог.
И попробуйте открыть где-нибудь ещё раз.
Блиннн. Похоже что вся проблема из-за того что не написал окончание (.tif).
Еще уточним Вы пишите:

Код: Выделить всё

r.out.gdal -f input=dif_jar@PERMANENT type=Float32 output=/home/lis/Dokumenty/dp/smaz/e/ИМЯ_ФАЙЛА.tif
флаг -f, но потом выскакивает ERROR все равно . А вот после флага -c , все в порядке. Какой тогда применить флаг? -c?

Ответить

Вернуться в «GRASS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость