Обработка Астрономических изображений из интернета

[Mercury127]

а, и у вас какой то нестандартный рандом...
вот у меня правильный, стандартный рандом:

G:\OLD_C\ICVCPP\BCC55\USER\RND>TYPE RND.C
#include <stdio.h>
#include <stdlib.h>

#define Rand(Range)    (rand()*(double)(Range)/RAND_MAX)

//То на 158 шаге получим слудующий адрес и таблицу

int i;
float x, y;

void main(){
        for( i= 0; i <160; i++){
                y=-90+Rand(180.0);
                x=Rand(360.0);
                //r=10;
                printf("%d %f %f\n", i, x, y);
        }// for
}// main

G:\OLD_C\ICVCPP\BCC55\USER\RND>..\..\Bin\bcc32.exe -I..\..\Include -L..\..\Lib RND.C
Borland C++ 5.5 for Win32 Copyright (c) 1993, 2000 Borland
rnd.c:
Turbo Incremental Link 5.00 Copyright (c) 1997, 2000 Borland
Warning: Unable to load DLL LNKDFM50

G:\OLD_C\ICVCPP\BCC55\USER\RND>RND.EXE > OUT

G:\OLD_C\ICVCPP\BCC55\USER\RND>TYPE OUT

0 120.655540 -89.285866
1 128.060547 -84.012268
2 193.310349 -50.903957
3 252.121948 -54.760277
4 98.923920 80.985443
5 39.233376 -10.028077
...
150 145.672165 -1.853999
151 232.752472 -73.487045
152 91.969360 -27.161779
153 46.374706 81.419418
154 94.111755 16.032290
155 197.902771 52.024597
156 52.131718 -5.303812
157 309.977722 23.486740
158 326.325867 27.123325
159 138.047424 83.545334

я даже поиск по файлу OUT сделал, нету нигде ваших 285.609314 -0.574053

[Finder BMap]

Спасибо за инфо, возможно это разница Borland/VS, немного совпадает только первое число в последовательности, я добавлю инфо в текст.

[Finder BMap]

Изображения  Неба от Телескопа Hubble

Теперь можно считывать файл изображений в формате FITS. Для этого используется следующий адрес, смысл суффикса "_c0m" можно найти в информации об ID, остальное ясно из контекста.

http://hla.stsci.edu/cgi-bin/getdata.cgi?config=ops&download=1&dataset=u66b0401r&filename=u66b0401r_c0m.fits


Полученный файл включает 4 изображения 800х800 в 32 bit float формате, соответсвующие 4 CCD инструмента, 3 из которых так скажем крупнозернистые, одно мелкозернистое (если я не ошибаюсь). Если сделать линейное преобразование в jpg Min…Max -> 0…255 то получиться следующее изображение(CCD #1). Синий цвет(линейный LUT с Blue Min =30) не является цветом обьектов, используется для уменьшения контрастности и усталости глаз.

[Finder BMap]

Изображения  Неба от Телескопа Hubble

Следующий шаг это очистка изображения от "артефактов". В данном случае это "нерабочие" зоны CCD удаляемые верхне-левым офсетом в 100 пикселов. Кроме то-го в FITS файлах могут встречаться пропущенные или пустые точки, помечаемые как NaN. Эти точки можно заменить 0 или что более корректно уровнем background, который находиться в параметрах данных. Проверка NaN в С++ может быть осуществлена одним из самых "логичных" в этом языке способов:

float Val ; if (Val==Val) then хорошая точка; else плохая точка;

что сравнимо с непоставлением "; " перед "else" в Pascal.
Можно также минимизировать данные этим уровнем, полагая что опуститься ниже сигнал уже не может. Если сделать какой-нибудь статистический Autoscale в jpeg, то может получиться примерно такое изображение. Данные готовы к собственно processing.

[Finder BMap]

Изображения  Неба от Телескопа Hubble

Некоторая дополнительная информация. Файл FITS содержит параметры данных, большинство из которых связаны с обстоятельствами наблюдения. Один из них background упомянут в предыдущем тексте. Извлечение параметров можно осуществить следующим образом
         
ffgkyd(infptr, "BACKGRND", &Val, NULL, &status);

Итерирование параметров:

fits_read_record(infptr, i, card, &status);

Полный список на обычных сайтах.
Кроме того включена информация об использовании изображений

"Based on observations made with the NASA/ESA Hubble Space Telescope,  and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). One copy of each paper resulting from data obtained from the HLA should be sent to the STScI."

Существуют также функции для координатных преобразований

ffgics, ffxypx, ffwldp etc.

[xd]

Горшочек, не вари!

[Finder BMap]

Изображения  Неба от Телескопа Hubble

2 немножко связанных с сюжетом вопроса для личной любознательности к знающим(извините за наивность если что)

1 Использовался ли Телескоп Hubble для научных наблюдений Земли?
2 Была прочитана информация, что американское военное ведомство передало NASA спутник-шпион с похожим телескопом. Правда ли это и если да использовался(зуется) ли он?

[xd]

1. Нет
2. С точностью до критериев похожести - нет.

[Finder BMap]

2. С точностью до критериев похожести - нет.

Спасибо, я понял Ваш ответ как передан но не такой. В инфо было написано вроде бы о размере зеркала. Если передан, Использовался(зуется) ли этот телескоп?

[xd]

Он - это кто?
И без конспирологии пожалуйста.

[Mercury127]

НАСА были то ли подарены, то ли проданы по дешевке два 2,4 м зеркала. одно использовали для Хабла, второе вроде собираются использовать в каком то новом спутнике хз когда.

для наблюдений Земли эти телескопы непригодны по то же простой причине, по которой ваши глаза непригодны для наблюдения Солнца в бинокль - ослепнут-с, понимаешь...

[Finder BMap]

Изображения  Неба от Обсерватории ESO

Изображения от Обсерватории ESO, затерянной в пустыне под южным солнцем, где как известно немного жарко но несказанно приятно, получаються сходным способом. Первый запрос теперь выглядит так

"http://archive.eso.org/wdb/wdb/eso/eso_archive_main/query?wdbo=csv%2fdisplay&max_rows_returned=200&instrument=&target=&resolver=simbad&ra=RARA&dec=DECDEC&box=00%20RADRAD%2000&degrees_or_hours=degrees&tab_target_coord=on&format=DecimDeg&wdb_input_file=&night=&stime=&starttime=12&etime=&endtime=12&tab_prog_id=on&prog_id=&gto=&pi_coi=&obs_mode=&title=&image[]=FORS2%3b%27IMA%25%27&image[]=HAWKI%3b%27IMA%25%27&image[]=VIMOS%3b%27IMA%25%27&tab_dp_cat=on&dp_cat=SCIENCE&tab_dp_type=on&dp_type=&dp_type_user=&tab_dp_tech=on&dp_tech=&dp_tech_user=&tab_dp_id=on&dp_id=&origfile=&tab_rel_date=on&rel_date=&obs_name=&ob_id=&tab_tpl_start=on&tpl_start=&tab_tpl_id=on&tpl_id=&tab_exptime=on&exptime=&tab_filter_path=on&filter_path=&gris_path=&grat_path=&slit_path=&tab_instrument=on&add=((ins_id%20like%20%27FORS2%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27))%20or%20(ins_id%20like%20%27HAWKI%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27))%20or%20(ins_id%20like%20%27VIMOS%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27)))&tab_tel_airm_start=on&tab_stat_instrument=on&tab_ambient=on&tab_stat_exptime=on&tab_HDR=on&tab_mjd_obs=on&aladin_colour=aladin_instrument&tab_stat_plot=on&order=&"

Где RARA, DECDEC, RADRAD заменяються на соответствующие величины.

Остальные параметры если не ясны из контекста, представлены на
http://archive.eso.org/eso/eso_archive_main.html

[Finder BMap]

Изображения  Неба от Обсерватории ESO

Если мы Демократически, без изменений применим туже случайную проседуру, то на 2059 шаге получим слудующий адрес и таблицу

http://archive.eso.org/wdb/wdb/eso/eso_archive_main/query?wdbo=csv%2fdisplay&max_rows_returned=200&instrument=&target=&resolver=simbad&ra=189.048&dec=25.9532&box=00%2010%2000&degrees_or_hours=degrees&tab_target_coord=on&format=DecimDeg&wdb_input_file=&night=&stime=&starttime=12&etime=&endtime=12&tab_prog_id=on&prog_id=&gto=&pi_coi=&obs_mode=&title=&image[]=FORS2%3b%27IMA%25%27&image[]=HAWKI%3b%27IMA%25%27&image[]=VIMOS%3b%27IMA%25%27&tab_dp_cat=on&dp_cat=SCIENCE&tab_dp_type=on&dp_type=&dp_type_user=&tab_dp_tech=on&dp_tech=&dp_tech_user=&tab_dp_id=on&dp_id=&origfile=&tab_rel_date=on&rel_date=&obs_name=&ob_id=&tab_tpl_start=on&tpl_start=&tab_tpl_id=on&tpl_id=&tab_exptime=on&exptime=&tab_filter_path=on&filter_path=&gris_path=&grat_path=&slit_path=&tab_instrument=on&add=((ins_id%20like%20%27FORS2%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27))%20or%20(ins_id%20like%20%27HAWKI%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27))%20or%20(ins_id%20like%20%27VIMOS%25%27%20AND%20(dp_tech%20like%20%27IMA%25%27)))&tab_tel_airm_start=on&tab_stat_instrument=on&tab_ambient=on&tab_stat_exptime=on&tab_HDR=on&tab_mjd_obs=on&aladin_colour=aladin_instrument&tab_stat_plot=on&order=&

OBJECT,RA,DEC,Program_ID,Instrument,Category,Type,Mode,Dataset ID,Release_Date,TPL ID,TPL START,Exptime,Exposure,Filter,MJD-OBS,Airmass,DIMM Seeing at Start
FORS OBJECT,189.087042,25.987556,60.A-9203(E),FORS2,SCIENCE,OBJECT,IMAGE,FORS2.2000-02-01T08:33:37.817,Jan 24 2001,FORS2_img_obs_exp,,   180.000,180.000,B_BESS,51575.356688,1.573,0.61
FORS OBJECT,189.087042,25.987556,60.A-9203(E),FORS2,SCIENCE,OBJECT,IMAGE,FORS2.2000-02-01T08:37:40.636,Jan 24 2001,FORS2_img_obs_exp,,    89.991,89.991,V_BESS,51575.359498,1.573,0.62
FORS OBJECT,189.087042,25.987556,60.A-9203(E),FORS2,SCIENCE,OBJECT,IMAGE,FORS2.2000-02-01T08:40:13.157,Jan 24 2001,FORS2_img_obs_exp,,    59.993,59.993,R_SPECIAL,51575.361263,1.574,0.53
...

Теперь используеться колонка 0 для имени и 8 для файла. Имя файла получается добавлением префикса "https://dataportal.eso.org/dataportal_new/file/"  :

https://dataportal.eso.org/dataportal_new/file/FORS2.2000-02-01T08:33:37.817


Возвращаемый файл ("FORS2.2000-02-01T08_33_37.817.fits.Z" в данном случае) может быть празднично упакован сначала компрессором FITS в .fz(иногда мнимым) затем UNIX в Z. Рекомендуется во всех сходных действиях соблюдать размер букв, поскольку если Вы звучите не в том регистре щепетильный сервер может не отдать файл. Распаковку можно сделать командной линией 7z и funpack.c для FTIS. Тем же или иным типом запроса можно добраться и до обработанных данных.

Результат может выглядеть следующим образом :

[xd]

Комментарий модератора

Хватит издеваться над фотографиями.
Вашу бы энергию да в мирное русло.

Тема закрыта