Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  377 1606 (N12E64 L=268)   493     96        3 ( 354/1584- 408/1624)  18.4 E67-E60 N09-N12
      2  236 1285 (S10E61 L=271)   723    104        1 ( 211/1228- 268/1333)  81.9 E62-E59 S13-S06
      3  190 1094 (S22E60 L=272)   774    167        0 ( 151/1041- 227/1146)  54.5 E63-E56 S25-S18
      4 1607 1050 (N12W33 L=  5)   300     78        3 (1591/1006-1621/1089) -53.1 W31-W34 N09-N14
      5  529  942 (S24E28 L=304)   530    114        0 ( 477/ 908- 591/ 977)  -6.3 E31-E22 S24-S22
      6 1743  720 (S00W55 L= 27)   452     47        1 (1720/ 697-1769/ 747)  90.0 W53-W57 S01-N01
      7 1156  475 (S33W27 L=359)   348     82        0 (1128/ 449-1191/ 513)  -9.5 W23-W29 S33-S31
      8 1357  386 (S31W44 L= 16)   938    189        3 (1283/ 356-1428/ 425)  14.0 W38-W50 S32-S28
  

    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