Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  388 1522 (N21E56 L=229)   601     91        1 ( 364/1480- 419/1562) -68.2 E57-E53 N18-N23
      2  284 1338 (N08E59 L=226)   775     89        1 ( 258/1298- 311/1382)  68.2 E60-E56 N05-N10
      3  926 1315 (N16E11 L=274)   486    189        0 ( 842/1268- 978/1349)  12.5 E15-E06 N12-N18
      4  306 1252 (N03E54 L=231)   560    113        0 ( 285/1195- 329/1301) -39.8 E57-E50 N00-N05
      5  749 1058 (S02E17 L=268)  1060    309        1 ( 685/ 933- 835/1193) -63.4 E21-E10 S08-N06
      6  514  972 (S11E32 L=253)   296    134        0 ( 492/ 916- 565/1001)  56.3 E33-E28 S14-S08
      7 1261  720 (S15W19 L=304)   202     46        0 (1235/ 704-1290/ 737)  45.0 W17-W21 S16-S13
      8 1633  440 (S27W58 L=342)   628    130        2 (1605/ 420-1663/ 467) -14.0 W53-W61 S28-S25
  

    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