Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1376 1461 (N33W18 L= 16)  2035    443        3 (1218/1329-1570/1565) -15.9 W04-W33 N26-N38
      2  928 1501 (N28E14 L=344)   424    129        1 ( 910/1446- 945/1564) -69.4 E16-E12 N24-N32
      3  388 1500 (N19E55 L=304)   495    127        0 ( 366/1459- 409/1534) -33.7 E57-E50 N17-N21
      4 1712 1459 (N39W48 L= 46)   410     70        0 (1686/1442-1738/1474)   0.0 W45-W50 N37-N39
      5  952 1278 (N14E08 L=350)   525    181        2 ( 875/1222-1028/1354) -26.6 E14-E02 N11-N17
      6  432 1291 (N07E43 L=315)   378     71        1 ( 414/1257- 450/1338)  90.0 E44-E41 N05-N10
      7 1830 1205 (N23W57 L= 55)  1110    210        0 (1771/1157-1881/1244) -13.0 W50-W63 N20-N25
      8 1343 1151 (N12W18 L= 16)   649    194        3 (1289/1061-1384/1239)  71.6 W15-W20 N06-N18
  

    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