private void LoadStars(SearchZone zone, double limitMag, List<IStar> starsFromThisZone) { List<LoadPosition> searchIndexes = m_Index.GetLoadPositions(zone); IntPtr buffer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(UCAC3Entry))); try { foreach (LoadPosition pos in searchIndexes) { string fileName; fileName = Path.Combine(m_CatalogLocation, string.Format("z{0}", pos.ZoneId.ToString("000"))); long positionFrom = (pos.FromRecordId - 1) * UCAC3Entry.Size; uint firstStarNoToRead = pos.FirstStarNoInBin + pos.FromRecordId; uint numRecords = pos.ToRecordId - pos.FromRecordId; using (FileStream str = new FileStream(fileName, FileMode.Open, FileAccess.Read, FileShare.Read)) { using (BinaryReader rdr = new BinaryReader(str)) { rdr.BaseStream.Position = positionFrom; for (int i = 0; i <= numRecords; i++, firstStarNoToRead++) { UCAC3Entry entry = new UCAC3Entry(); byte[] rawData = rdr.ReadBytes(UCAC3Entry.Size); Marshal.Copy(rawData, 0, buffer, Marshal.SizeOf(entry)); entry = (UCAC3Entry)Marshal.PtrToStructure(buffer, typeof(UCAC3Entry)); entry.InitUCAC3Entry((ushort)pos.ZoneId, firstStarNoToRead); if (entry.Mag > limitMag) continue; if (entry.RAJ2000 < zone.RAFrom) continue; if (entry.RAJ2000 > zone.RATo) continue; if (entry.DEJ2000 < zone.DEFrom) continue; if (entry.DEJ2000 > zone.DETo) continue; starsFromThisZone.Add(entry); } } } } } finally { Marshal.FreeHGlobal(buffer); } }
private void LoadStars(SearchZone zone, double limitMag, List<IStar> starsFromThisZone) { List<LoadPosition> loadPositions = m_Index.GetLoadPositions(zone); foreach (LoadPosition loadPos in loadPositions) { LoadStarsBrighterThan(zone, loadPos, limitMag, starsFromThisZone); } }
private void LoadStarsBrighterThan(SearchZone zone, LoadPosition loadPos, double limitMag, List<IStar> starsFromThisZone) { IntPtr buffer = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(NOMADEntry))); try { using (FileStream fs = new FileStream(m_Index.GetCatalogFileForZone(loadPos.ZoneId), FileMode.Open, FileAccess.Read)) using (BinaryReader rdr = new BinaryReader(fs)) { #if !PRODUCTION Trace.Assert(fs.Length >= (loadPos.FromRecordId + 1) * 88); #endif uint currStarInBin = loadPos.FromRecordId - 1; fs.Position = 88 * loadPos.FromRecordId; for (int i = 0; i < loadPos.ToRecordId - loadPos.FromRecordId; i++) { currStarInBin++; #if !PRODUCTION Trace.Assert(fs.Position % 88 == 0); #endif fs.Seek(44, SeekOrigin.Current); double mag = rdr.ReadInt32() / 1000.0; /* V Mag */ int byteDiff = 0; if (mag == 30000.0) { mag = rdr.ReadInt32() / 1000.0; /* R Mag */ byteDiff = 4; } if (mag > limitMag) { fs.Seek(40 - byteDiff, SeekOrigin.Current); continue; } fs.Seek(-(48 + byteDiff), SeekOrigin.Current); NOMADEntry entry = new NOMADEntry(); byte[] rawData = rdr.ReadBytes(NOMADEntry.Size); Marshal.Copy(rawData, 0, buffer, Marshal.SizeOf(entry)); entry = (NOMADEntry)Marshal.PtrToStructure(buffer, typeof(NOMADEntry)); if (entry.RAJ2000 >= zone.RAFrom && entry.RAJ2000 <= zone.RATo) { entry.InitNOMADEntry(loadPos.FirstStarNoInBin + currStarInBin, (uint)loadPos.ZoneId, currStarInBin); starsFromThisZone.Add(entry); } } } } finally { Marshal.FreeHGlobal(buffer); } }
public List<LoadPosition> GetLoadPositions(SearchZone zone) { List<LoadPosition> loadPositions = new List<LoadPosition>(); int deltaFrom = zone.DEFrom > 0 ? 1 : 0; int deltaTo = zone.DETo > 0 ? 1 : 0; int fromZone = Math.Max(1, (int)(zone.DEFrom * 5) + 90 * 5 + deltaFrom); int toZone = Math.Min(900, (int)(zone.DETo * 5) + 90 * 5 + deltaTo); for (int zoneId = fromZone; zoneId <= toZone; zoneId++) { LoadPosition loadPos = new LoadPosition(); loadPos.ZoneId = zoneId; loadPos.FirstStarNoInBin = 0; int raFrom = Math.Max(0, (int)Math.Floor(zone.RAFrom)); int raTo = Math.Min(360, (int)Math.Ceiling(zone.RATo) + 1); loadPos.FromRecordId = (uint)ZoneIndex[zoneId - 1].RAStartPositions[raFrom]; if (raTo == 360) loadPos.ToRecordId = (uint)ZoneIndex[zoneId - 1].TotalStarsInZone; else loadPos.ToRecordId = (uint)ZoneIndex[zoneId - 1].RAStartPositions[raTo]; loadPositions.Add(loadPos); } return loadPositions; }
internal List<LoadPosition> GetLoadPositions(SearchZone zone) { List<LoadPosition> loadPositions = new List<LoadPosition>(); int signDEFrom = Math.Sign(zone.DEFrom); int fromZone, toZone; bool isSouth = signDEFrom < 0; fromZone = (int)Math.Floor(((zone.DEFrom + 90) * 4)); toZone = (int)Math.Ceiling(((zone.DETo + 90) * 4)); for (int zoneId = fromZone; zoneId < toZone; zoneId++) { LoadPosition loadPos = new LoadPosition(); if (zoneId == 359) { loadPos.ZoneId = 0; loadPos.SubZoneId = 'a'; loadPos.Hemisphere = 's'; } else if (zoneId == 360) { loadPos.ZoneId = 0; loadPos.SubZoneId = 'a'; loadPos.Hemisphere = 'n'; } else { if (isSouth) { loadPos.ZoneId = 89 - (zoneId / 4); loadPos.SubZoneId = "dcba"[zoneId % 4]; loadPos.Hemisphere = 's'; } else { loadPos.ZoneId = (zoneId / 4) - 90; loadPos.SubZoneId = "abcd"[zoneId % 4]; loadPos.Hemisphere = 'n'; } } int raFrom = Math.Max(0, (int)Math.Floor(zone.RAFrom)); int raTo = Math.Min(359, (int)Math.Ceiling(zone.RATo)); loadPos.FromRecordId = (uint)ZoneIndex[zoneId].RAIndex[raFrom]; loadPos.ToRecordId = (uint)ZoneIndex[zoneId].RAIndex[raTo]; loadPositions.Add(loadPos); } return loadPositions; }
private void LoadStars(SearchZone zone, double limitMag, List<IStar> starsFromThisZone) { List<LoadPosition> searchIndexes = m_Index.GetLoadPositions(zone); foreach (LoadPosition pos in searchIndexes) { string fileName; fileName = Path.Combine(m_CatalogLocation, string.Format("{0}{1}{2}.dat", pos.Hemisphere, pos.ZoneId.ToString("00"), pos.SubZoneId)); long positionFrom = (pos.FromRecordId + 1 /* for the header row */) * PPMXLEntry.Size; uint numRecords = pos.ToRecordId - pos.FromRecordId; using (FileStream fileStr = new FileStream(fileName, FileMode.Open, FileAccess.Read)) using (BinaryReader rdr = new BinaryReader(fileStr)) { rdr.BaseStream.Position = positionFrom; for (int i = 0; i <= numRecords; i++) { byte[] data = rdr.ReadBytes(PPMXLEntry.Size); PPMXLEntry entry = new PPMXLEntry(Encoding.ASCII.GetString(data)); if (entry.Mag > limitMag) continue; if (entry.RAJ2000 < zone.RAFrom) continue; if (entry.RAJ2000 > zone.RATo) continue; if (entry.DEJ2000 < zone.DEFrom) continue; if (entry.DEJ2000 > zone.DETo) continue; starsFromThisZone.Add(entry); } } } }