เบื้องหลังการถ่ายทำ steetview ของกูเกิ้ล

มาดูกันชัดๆว่ารถเก็บข้อมูล google streetview ทำงานอย่างไร พนักงานที่กูเกิ้ลประกาศรับสมัครงาน ที่ว่าต้องสามารถปีนขึ้นหลังคารถและติดตั้งเครื่องมือประจำวันได้ทำงานหนักอย่างไร

Google Street View บุกประเทศไทย

แล้วข่าวที่คาดเดาก็เปิดเผยออกมาแล้ว เมื่อกูเกิ้ลแถลงข่าวร่วมพันธมิตรคือ การท่องเที่ยวแห่งประเทศไทย (ททท.) เริ่มโครงการสตรีทวิวที่ทั่วประเทศไทย "เพื่อสนับสนุนแคมเปญใหม่ในการรณรงค์การท่องเที่ยวไทย" ยินดีด้วยครับ
-------------------------------------------------------
---------------------กำหนดการแถลงข่าว----------------
--------“ส่งเสริมการท่องเที่ยวไทย ด้วย Google แผนที่กับสตรีทวิว”-----------
-------เพื่อประกาศเริ่มโครงการ Google แผนที่กับสตรีทวิว ทั่วประเทศ-----------
------------วันพฤหัสบดีที่ ๑ กันยายน ๒๕๕๔ ------------------------
----------ณ ร้านกาแฟดอยตุง พระที่นั่งอนันตสมาคม-------------------
------------------------------------------------------

ฝูงของรถกูเกิ้ลสตรีทวิว นับได้ 15 คัน

GPU กับงานประมวลผลภูมิสารสนเทศ

ขอบคุณดร.ฮอนด้า(AIT) ดร.ภุชงค์ (KU) และ ศาสตราจารย์ อาโอกิ (Tokyo Univ.) จัดฝึกอบรม
Joint International Workshop on HPC for Natural Disaster Simulation and GPU computing
Dr. Putchong Uthayopas, Head of computer engineering, Kasetsart University
Prof. Aoki Takayuki, GSIC, Tokyo Institute of Technology
Prof. Honda Kiyoshi, IDEAS Chubu University

27 June 2011: Presentation Session
28 June 2011: GPU Tutorial Session

คู่มือการประมวลผลภูมิสารสนเทศด้วย FOSS4G

คู่มือการประมวลผลภูมิสารสนเทศด้วย FOSS4Gล่าสุดครับ

รหัสมาตรฐานการแปลงพื้นหลักฐานและการฉายแผนที่

ในบ้านเรา รหัสมาตรฐานการแปลงพื้นหลักฐานและการฉายแผนที่ ที่ต้องใช้บ่อย ผมทบทวนให้ดูดังนี้
มาตรฐานรหัสพื้นหลักฐานแผนที่ (map datum) และการฉายแผนที่ (map projection) ตามมาตรฐาน European Petroleum Survey Group (EPSG) ปัจจุบันชื่อว่าThe OGP Surveying and Positioning Committee มาตรฐานนี้เป็นที่ยอมรับมาตรฐานานาชาติ ISO และมาตรฐานวิชาชีพ OGC

พื้นหลักฐาน	ระบบพิกัดและการฉายแผนที่	รหัสมาตรฐาน
WGS84	Geographic Coord.Sys.	EPSG:4326
WGS84	UTM Zone 47N	EPSG:32647
WGS84	UTM Zone 48N	EPSG:32648
IndianDatum1975	Geographic Coord.Sys.	EPSG:4240
IndianDatum1975	UTM Zone 47N	EPSG:24047
IndianDatum1975	UTM Zone 48N	EPSG:24048

*** ข้อสังเกตและข้อควรระวัง พืนหลักฐาน IndianDatum1975 กำหนดขนาดของทรงรีโลกและพารามิเตอร์ต่างๆชัดเจน แต่!!!ไม่รวมถึงพารามิเตอร์การแปลงพื้นหลักฐาน ผู้ใช้จะต้องเลือกใช้เอง , สรุปสั้นๆ ไม่มี "ตัวเลข" มาตรฐานพารามิเตอร์การแปลงพื้นหลักฐานะครับ แต่มีชื่อฟิลด์ให้ใช้ครับ เช่น มักจะเรียกว่า TOWGS84 ครับ ***

กรณีพื้นที่ศึกษาของท่านอยู่ในพื้นที่อื่นๆนอกเหนือจากประเทศไทย ท่านสามารถตรวจสอบรหัส CRS ได้ที่ http://www.epsg.org/ or http://spatialreference.org/?

การแปลงพื้นหลักฐานและการฉายแผนที่ด้วย FOSS4G

มาตรฐานการแปลงพื้นหลักฐานและการฉายแผนที่ เป็นเรื่องสำคัญและยังสร้างความสับสนกับวงการแผนที่ไทย โดยเฉพาะซอฟต์แวร์ที่โดยธรรมชาติไม่ได้บรรจุพารามิเตอร์การแปลงที่ประเทศต้องใช้ΔX= -206 m., ΔY=-837 m., ΔZ=-295 m.
โชดดีที่ซอฟต์แวร์ไลบรารี่ proj4 มีความสดวกนี้รองรับ ซึ่งมีการเรียกใช้ในซอฟต์แวร์ gdalwarp สำหรับการแปลงแผนที่ภาพ และ ogr2ogr สำหรับการแปลงเวกเตอร์

ก่อนทำการแปลงพื้นหลักฐานและการฉาย ผู้ใช้ต้องเตรียมไฟล์สั้นที่เป็นนิยามพื้นหลักฐานและการฉายดังตารางนี้

----------------------------------------------------------------------
ชื่อไฟล์ - พื้นหลักฐาน - นิยามของพื้นหลักฐานและการฉาย
----------------------------------------------------------------------
ThaiInd1975.pj4
พื้นหลักฐานอินเดียน 1975และระบบพิกัดยีออเดติกส์
+proj=lonlat +ellps=evrst30 +towgs84=206,837,295
----------------------------------------------------------------------
ThaiInd1975_utm47.pj4
พื้นหลักฐานอินเดียน 1975และการฉายระบบพิกัดยูทีเอ็มโซน47
+proj=tmerc +ellps=evrst30 +lon_0=99 +x_0=500000 +k=0.9996 +towgs84=206,837,295
----------------------------------------------------------------------
ThaiInd1975_utm48.pj4
พื้นหลักฐานอินเดียน 1975 และการฉายระบบพิกัดยูทีเอ็มโซน48
+proj=tmerc +ellps=evrst30 +lon_0=105 +x_0=500000 +k=0.9996 +towgs84=206,837,295
----------------------------------------------------------------------

+++ตัวอย่างสำหรับการฉายแผนที่ภาพแรสเตอร์ ด้วยยูทิลิตี้ gdalwarp ในการแปลงภาพจากพื้นหลักฐานอินเดียน 1975 ระบพิกัดยูทีเอ็มโซน 47 ให้เป็นแผนที่ภาพใหม่ที่ระบบพิกัดอ้างอิงบนพื้นหลักฐาน WGS84 และให้ระบบพิกัดเป็นค่าพิกัดยีออเดติกส์ (ระบบพิกัดภูมิศาสตร์) เพื่อเตรียมไปซ้อนทับบน Google Earth

gdalwarp -s_srs “ThaiInd75_utm47.pj4” -t_srs epsg:4326 ortho_id75utm47.tif ortho_geo_wgs84.tif

+++ มีข้อมูลเวกเตอร์ในระบบพื้นหลักฐานอินเดียน 1975 ฉายด้วยระบบพิกัดยูทีเอ็มโซน 47 ต้องการแปลงให้เป็นพื้นหลักฐาน WGS84 และให้ระบบค่าพิกัดเป็นค่าพิกัดยีออเดติกส์ (ภูมิศาสตร์) สำหรับนำไปใช้งานบน Google Earth

ogr2ogr –f “KML” -t_srs EPSG:4326 -s_srs “ThaiInd1975_utm47.pj4” province_geo_wgs84.kml province_id75utm47.shp

PL/PgSQL สำหรับสร้าง 'ช่วงเวลาของวัน'

นานมาแล้ว ได้เขียน PL/PgSQL สำหรับสร้าง 'ช่วงเวลาของวัน' สำหรับการวิเคราะห์ช่วงเวลาของการเกิดเหตุการณ์ โดยผู้ใช้กำหนดว่า
-6โมงเช้าถึง 7โมงเช้าเรียกว่า สล็อต1
-7โมงเช้าถึง 8โมงเช้าเรียกว่า สล็อต2
.....
และเพื่อใช้ในการ query ก็เลยเขียน PL ไว้ แต่กลัวหาย และอาจเป็นประโยชน์ ก็นำมาฝากไว้ที่นี่

ตัวอย่างการใช้งาน
FUNCTION "PS_Hour_6Slot24"(datetime timestamp with time zone, slot_number boolean)



==================== โค้ดตามนี้ครับ ===========================
-- Function: "PS_Hour_6Slot24"(timestamp with time zone, boolean)

-- DROP FUNCTION "PS_Hour_6Slot24"(timestamp with time zone, boolean);

CREATE OR REPLACE FUNCTION "PS_Hour_6Slot24"(datetime timestamp with time zone, slot_number boolean)
RETURNS text AS
$BODY$
DECLARE
hour_int integer;
hour_int_1 integer;
hour_int_2 integer;
slot_text text;
BEGIN
IF slot_number = TRUE THEN
hour_int := (to_char( datetime, 'HH24' )::integer);
hour_int := hour_int - 5;
IF hour_int < 0 THEN
hour_int := hour_int + 24;
ELSEIF hour_int = 0 THEN
hour_int := 24;
END IF;
RETURN hour_int;
ELSE
hour_int_1 := (to_char( datetime, 'HH24' )::integer);
hour_int_2 := hour_int_1 + 1;
IF hour_int_2 = 24 THEN
hour_int_2 := 0;
END IF;
slot_text := (to_char(hour_int_1, 'FM00')) || '-' || (to_char(hour_int_2, 'FM00' )) ;
RETURN slot_text;
END IF;


END;$BODY$
LANGUAGE plpgsql VOLATILE
COST 100;
ALTER FUNCTION "PS_Hour_6Slot24"(timestamp with time zone, boolean) OWNER TO postgres;
COMMENT ON FUNCTION "PS_Hour_6Slot24"(timestamp with time zone, boolean) IS 'return either hour-slot 1 to 24 starting from 6 o''clock or hour-slot text';
==================================================================================