Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  295 1194 (S13E52 L=357)  1593    316        3 ( 241/1097- 366/1293)  76.0 E57-E46 S18-S06
      2  646 1033 (S15E22 L= 27)   444    180        0 ( 593/ 957- 688/1121) -37.9 E27-E18 S19-S11
      3 1713 1068 (N14W41 L= 91)   334     86        2 (1687/1034-1742/1117) -78.7 W39-W43 N11-N16
      4 1326  978 (S02W18 L= 67)   957    245        0 (1223/ 941-1438/1054)   0.0 W10-W24 S04-N00
      5 1852 1003 (N15W56 L=105)   457    150        0 (1822/ 981-1882/1037)  38.7 W52-W58 N13-N18
      6 1757  838 (N03W52 L=101)   955    168        1 (1724/ 761-1786/ 913) -60.3 W48-W56 N00-N07
      7  406  844 (S31E37 L= 13)   816    215        1 ( 312/ 816- 493/ 867)  11.3 E44-E29 S34-S29
      8 1156  727 (S20W15 L= 65)   236     85        0 (1109/ 713-1202/ 741)  18.4 W12-W18 S22-S19
      9 1660  563 (S15W55 L=104)  1242    194        3 (1581/ 500-1721/ 651)  70.0 W50-W58 S19-S08
     10 1141  571 (S30W20 L= 69)   574    167        0 (1065/ 543-1234/ 613)  45.0 W15-W24 S33-S24
     11 1707  513 (S15W63 L=112)  1595    247        3 (1643/ 472-1765/ 559)  19.7 W57-W72 S19-S12
  

    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