Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1627 1360 (N33W31 L=124)   522     75        2 (1606/1324-1647/1402) -45.0 W28-W33 N30-N35
      2 1808 1100 (N23W52 L=145)  2586    394        4 (1700/ 970-1907/1187)  38.0 W43-W66 N12-N30
      3  914 1112 (S02E08 L= 85)   229     55        0 ( 888/1087- 935/1136) -33.7 E10-E07 S03-S01
      4  402  772 (S35E37 L= 57)  2552    364        8 ( 315/ 667- 492/ 914)  49.6 E47-E30 S44-S24
      5 1757  651 (S04W62 L=155)   536     27        0 (1738/ 624-1774/ 680) -33.7 W58-W64 S05-S02
      6 1381  596 (S18W34 L=128)   674    180        2 (1309/ 566-1449/ 658)  56.3 W30-W36 S21-S12
      7 1160  602 (S24W21 L=114)   354     84        0 (1126/ 583-1203/ 626)   0.0 W17-W23 S25-S22
      8 1451  554 (S18W41 L=134)   591    130        0 (1397/ 514-1503/ 583)  60.3 W38-W44 S22-S15
      9 1023  522 (S33W16 L=109)   782    227        2 ( 946/ 463-1120/ 581)  63.4 W11-W19 S39-S27
  

    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