Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1150 1875 (N59W04 L=302)   504     70        0 (1114/1861-1184/1887)   0.0 W00-W07 N56-N60
      2  798 1857 (N53E33 L=265)   597     52        1 ( 775/1827- 823/1882) -74.1 E36-E29 N49-N56
      3  992 1862 (N55E13 L=285)   509     77        1 ( 955/1840-1034/1879)  14.0 E16-E08 N52-N57
      4 1129 1643 (N36W03 L=300)   272     61        0 (1105/1608-1145/1674) -68.2 W00-W04 N33-N38
      5  380 1623 (N31E59 L=239)   755     61        1 ( 352/1594- 401/1650)  56.3 E63-E55 N28-N33
      6  733 1565 (N27E24 L=274)   667    181        1 ( 675/1511- 814/1588)  24.0 E28-E18 N23-N29
      7 1405 1450 (N25W22 L=320)   283    105        0 (1356/1423-1450/1475) -15.9 W18-W25 N23-N25
      8  193  983 (S10E60 L=237)   420     58        0 ( 177/ 947- 211/1018)  63.4 E62-E58 S12-S08
      9  904  831 (S17E06 L=291)   278    100        0 ( 884/ 779- 931/ 870)  63.4 E07-E04 S20-S14
     10  628  743 (S24E24 L=273)   580    158        0 ( 567/ 702- 706/ 769)  20.0 E29-E18 S26-S21
     11 1010  481 (S39W03 L=301)   848    146        3 ( 964/ 440-1061/ 526) -15.9 W00-W07 S41-S35
     12 1564  375 (S40W56 L=354)   390     73        1 (1527/ 362-1606/ 391)  23.2 W51-W60 S41-S38
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y