PsXisfReader.ps1
|
class XisfFileStats { [decimal]$Exposure [string]$Filter [string]$Instrument [string]$Object [decimal]$Gain [decimal]$Offset [string]$ImageType [decimal]$CCDTemp [decimal]$SetTemp [decimal]$FocalLength [decimal]$FocalRatio [nullable[DateTime]]$ObsDate [nullable[DateTime]]$ObsDateMinus12hr [nullable[DateTime]]$LocalDate [decimal]$SSWeight [decimal]$Pedestal [string[]]$History [System.IO.FileInfo]$Path [decimal]$XPIXSZ [decimal]$YPIXSZ [string]$Geometry [decimal]$Rotator [bool]HasTokensInPath([string[]]$tokens){ $hasToken=$false foreach( $x in $tokens) { if($this.Path.FullName.ToLower().Contains($x.ToLower())){ $hasToken=$true break; } } return $hasToken; } [bool]IsIntegratedFile() { return [bool]($this.History) } } Add-Type -AssemblyName 'System.Xml.Linq' Function Read-XisfSignature([System.IO.BinaryReader]$reader) { $buffer=[char[]]::new(8) $reader.Read($buffer,0,8)>$null [System.String]::new($buffer) } Function Read-XisfHeader([System.IO.BinaryReader]$reader) { $length = $reader.ReadUInt32() $buffer=[char[]]::new(4) $reader.Read($buffer,0,4)>$null $header=[byte[]]::new($length) $reader.Read($header,0,$length)>$null $memoryStream=[System.IO.MemoryStream]::new($header) try { $xmlReader=[System.Xml.XmlReader]::Create($memoryStream) try { [xml][System.Xml.Linq.XDocument]::Load($xmlReader) } finally { $xmlReader.Dispose() } } finally { $memoryStream.Dispose() } } Function Get-XisfImageScale { [CmdletBinding()] param( [Parameter(ValueFromPipeline=$true,Mandatory=$true)]$Stats ) Write-Output (206.265*$Stats.XPIXSZ/$Stats.FocalLength) } Function Get-XisfHeader { [CmdletBinding()] [OutputType([xml])] param ( [Parameter(ValueFromPipeline=$true,Mandatory=$true)][System.IO.FileInfo]$Path ) process{ $stream = [System.IO.File]::OpenRead($Path.FullName) $reader=new-object System.IO.BinaryReader $stream try { Read-XisfSignature $reader > $null Read-XisfHeader $reader } finally { $reader.Dispose() $stream.Dispose() } } } Function Get-XisfFitsStats { [CmdletBinding()] [OutputType([XisfFileStats])] param ( [Parameter(ValueFromPipeline=$true,Mandatory=$true)][System.IO.FileInfo]$Path, [Parameter(ValueFromPipeline=$false)][Hashtable]$Cache, [Parameter(ValueFromPipeline=$false)][Hashtable]$ErrorCache ) process{ if($Cache -and $Cache.ContainsKey($Path.FullName)) { $result = $Cache[$Path.FullName] Write-Output ([XisfFileStats]@{ Exposure=$result.Exposure Filter=$result.Filter Instrument=$result.Instrument Object=$result.Object Gain=$result.Gain Offset=$result.Offset ImageType=$result.ImageType CCDTemp=$result.CCDTemp SetTemp=$result.SetTemp FocalLength=$result.FocalLength FocalRatio=$result.FocalRatio ObsDate=$result.ObsDate ObsDateMinus12hr=$result.ObsDateMinus12hr LocalDate=$result.LocalDate SSWeight=$result.SSWeight Pedestal=$result.Pedestal History=$result.History Path=$Path XPIXSZ=$result.XPIXSZ YPIXSZ=$result.YPIXSZ Geometry=$result.Geometry Rotator=$result.Rotator }) } elseif($ErrorCache -and $ErrorCache.ContainsKey($Path.FullName)){ Write-Verbose "Skipping file $($Path.FullName) due to previously encountered error." } else { $result=$null try { $header = Get-XisfHeader -Path $Path $geometry = $header.xisf.Image.geometry $fits = $header.xisf.Image.FITSKeyword $filter = $fits.Where({$_.Name -eq 'FILTER'}).value if($filter) { $filter = $filter.Replace('5nm','').Replace('3nm','').Replace("'","") } $obsDate = $fits.Where{$_.Name -eq 'DATE-OBS'}.value | Select-Object -First 1 $obsDateMinus12hr= $null if($obsDate){ $obsDateMinus12hr=$obsDate| foreach-object{ ([DateTime]($_.Trim("'"))).AddHours(-12).Date } } $results=@{ Exposure=$fits.Where{$_.Name -eq 'EXPOSURE'}.value Filter=$filter Instrument=$fits.Where{$_.Name -eq 'INSTRUME'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 Object=$fits.Where{$_.Name -eq 'OBJECT'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 Gain=$fits.Where{$_.Name -eq 'GAIN'}.value | select-object -First 1 Offset=$fits.Where{$_.Name -eq 'OFFSET'}.value | select-object -First 1 ImageType=$fits.Where{$_.Name -eq 'IMAGETYP'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 CCDTemp=$fits.Where{$_.Name -eq 'CCD-TEMP'}.value | select-object -First 1 SetTemp=$fits.Where{$_.Name -eq 'SET-TEMP'}.value | select-object -First 1 FocalLength=$fits.Where{$_.Name -eq 'FOCALLEN'}.value|%{if($_){$_.Trim("'").TrimEnd("mm")}} | select-object -First 1 FocalRatio=$fits.Where{$_.Name -eq 'FOCRATIO'}.value | select-object -First 1 ObsDate=$obsDate|%{if($_){$_.Trim("'")}} | select-object -First 1 ObsDateMinus12hr=$obsDateMinus12hr LocalDate=$fits.Where{$_.Name -eq 'DATE-LOC'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 SSWeight=$fits.Where{$_.Name -eq 'SSWEIGHT'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 Pedestal=$fits.Where{$_.Name -eq 'PEDESTAL'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 History=$fits.Where{$_.Name -eq 'HISTORY'}.comment XPIXSZ=$fits.Where{$_.Name -eq 'XPIXSZ'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 YPIXSZ=$fits.Where{$_.Name -eq 'YPIXSZ'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 Rotator=$fits.Where{$_.Name -eq 'ROTATOR'}.value|%{if($_){$_.Trim("'")}} | select-object -First 1 Path=$Path Geometry=$geometry } $result = [XisfFileStats] $results if($Cache) { $Cache.Add($Path.FullName,$result) } } catch [System.Xml.XmlException]{ Write-Warning ("An error occured reading the file "+($Path.FullName)) Write-Verbose $_.Exception.ToString() if($ErrorCache){ $ErrorCache.Add($Path.FullName,$_.Exception) } } catch { Write-Warning "An unexpected error occured processing file $($Path.FullName)" Write-Warning $_.Exception.ToString() if($ErrorCache){ $ErrorCache.Add($Path.FullName, $_.Exception) } throw } if($result){ Write-Output $result } } } } Function Get-XisfProperties { [CmdletBinding()] param ( [Parameter(ValueFromPipeline=$true,Mandatory=$true)][System.IO.FileInfo]$Path ) process{ $stream = [System.IO.File]::OpenRead($Path.FullName) $reader=new-object System.IO.BinaryReader $stream try { $header=Get-XisfHeader -Path $Path $header.Xisf.Image.Property } catch [System.Xml.XmlException]{ Write-Warning ("An error occured reading the file "+($Path.FullName)) Write-Verbose $_.Exception.ToString() if($ErrorCache){ $ErrorCache.Add($Path.FullName,$_.Exception) } } catch { Write-Warning "An unexpected error occured processing file $($Path.FullName)" Write-Warning $_.Exception.ToString() throw } finally { $reader.Dispose() $stream.Dispose() } } } Function Measure-ExposureTime { param( [Parameter(ValueFromPipeline=$true)][XisfFileStats[]]$Input, [Parameter()][Switch]$TotalSeconds, [Parameter()][Switch]$TotalMinutes, [Parameter()][Switch]$TotalHours, [Parameter()][Switch]$TotalDays, [Parameter()][Switch]$First, [Parameter()][Switch]$Last ) begin { $earliest=$null $latest=$null $seconds=0.0 $result=new-object psobject -Property @{ Count=0 } } process{ if($null -eq $earliest -or $earliest.ObsDate -gt $Input.ObsDate){ $earliest=$Input } if($null -eq $latest -or $earliest.ObsDate -lt $Input.ObsDate){ $latest=$Input } $seconds+=$Input.Exposure $result.Count+=1 } end { $total = [TimeSpan]::FromSeconds($seconds) if($TotalSeconds){ add-member -MemberType NoteProperty -InputObject $result -Name "TotalSeconds" -value $total.TotalSeconds } if($TotalMinutes){ add-member -MemberType NoteProperty -InputObject $result -Name "TotalMinutes" -value $total.TotalMinutes } if($TotalHours){ add-member -MemberType NoteProperty -InputObject $result -Name "TotalHours" -value $total.TotalHours } if($TotalDays){ add-member -MemberType NoteProperty -InputObject $result -Name "TotalDays" -value $total.TotalDays } if($First){ add-member -MemberType NoteProperty -InputObject $result -Name "First" -value $earliest } if($Last){ add-member -MemberType NoteProperty -InputObject $result -Name "Last" -value $latest } add-member -MemberType NoteProperty -InputObject $result -PassThru -Name "Total" -Value $total } } Function Get-MoonPercentIlluminated { [CmdletBinding()] param( [Parameter(Mandatory,ValueFromPipeline)][DateTime]$Date ) begin{ $julianConst=2415018.5 $knownDateOfNewMoon= [DateTime]::SpecifyKind( [DateTime]"1920-01-21 05:25:00", [DateTimeKind]::Utc).ToOADate() $moonDaysPerCycle=29.5306 } process{ $julianDate=$Date.ToUniversalTime().ToOADate()+$julianConst $daysSinceLastNewMoon=$knownDateOfNewMoon+$julianConst $newMoons = ($julianDate - $daysSinceLastNewMoon) / $moonDaysPerCycle; $daysIntoCycle = ($newMoons - [Math]::Truncate($newMoons))*$moonDaysPerCycle $x=($daysIntoCycle/$moonDaysPerCycle)*[Math]::PI [Math]::Round( [Math]::Sin($x),2) } } |