Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1158 1662 (N38W31 L=157)   716    129        2 (1137/1600-1180/1712)  36.9 W27-W35 N33-N40
      2 1533 1666 (N25W56 L=182)  1596    236        2 (1444/1625-1609/1694) -36.0 W49-W62 N21-N29
      3  653 1308 (N31E15 L=111)   800    359        2 ( 504/1284- 825/1333) -24.4 E25-E03 N25-N36
      4 1632 1129 (S06W39 L=165)  2235    650        7 (1554/ 845-1749/1290)  49.4 W27-W51 S20-N01
      5 1782 1207 (S07W54 L=180)   350     93        1 (1753/1183-1807/1233)   0.0 W50-W56 S07-S04
      6  969 1157 (N14W00 L=127)   363     64        1 ( 939/1126-1002/1187) -76.0 E01-W01 N12-N16
      7 1810 1056 (S16W51 L=178)   748    142        3 (1752/1020-1849/1097)  -5.2 W45-W56 S17-S14
      8 1274  767 (S16W07 L=133)   447    136        1 (1198/ 730-1332/ 797) -49.4 W03-W09 S19-S12
      9  809  485 (S20E28 L= 98)   288     90        1 ( 778/ 462- 831/ 519)  21.8 E31-E26 S21-S17
     10  959  403 (S29E23 L=103)   878    248        1 ( 928/ 310- 986/ 499)  61.7 E27-E19 S35-S22
     11  482  362 (S20E58 L= 68)   541    107        1 ( 463/ 333- 515/ 389)   8.1 E61-E54 S21-S18
     12 1140  252 (S47E20 L=106)  1470    233        1 (1058/ 206-1188/ 318)   0.0 E26-E14 S51-S41
  

    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