Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  563 1745 (N25E58 L=217)  1529    181        3 ( 508/1661- 607/1818) -49.8 E67-E54 N17-N30
      2  758 1762 (N30E43 L=231)   758    116        3 ( 728/1707- 785/1813) -35.0 E49-E39 N26-N33
      3  582 1675 (N20E50 L=225)  1720    290        4 ( 504/1595- 671/1737)  77.0 E53-E44 N13-N26
      4 1149 1631 (N32E11 L=264)   582    153        1 (1128/1557-1175/1704) -76.0 E13-E08 N26-N38
      5 1575 1476 (N36W23 L=298)   647    106        1 (1535/1433-1596/1522)  60.3 W20-W24 N33-N40
      6  397 1394 (S01E50 L=225)   614     88        1 ( 365/1357- 426/1444) -23.2 E54-E46 S03-N01
      7  245 1228 (S13E59 L=216)  1835    280        4 ( 183/1145- 308/1351)  72.9 E65-E55 S18-S05
      8  187 1064 (S24E63 L=212)   423     61        2 ( 155/1050- 224/1079)  11.3 E67-E57 S24-S22
      9 1250  890 (S08W16 L=291)   228     72        1 (1219/ 862-1274/ 924)  78.7 W15-W17 S10-S05
     10 1204  529 (S29W27 L=302)   934    231        1 (1104/ 505-1315/ 552)  28.3 W20-W33 S31-S24
     11 1502  517 (S20W47 L=322)   336    112        1 (1459/ 486-1533/ 546)   0.0 W42-W50 S21-S19
     12 1705  450 (S16W71 L=346)   670     89        1 (1653/ 425-1742/ 480)  39.8 W66-W75 S18-S13
  

    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