Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  652 1342 (N15E28 L=331)   583    263        1 ( 596/1223- 688/1432) -47.7 E34-E24 N08-N20
      2  337 1341 (N11E53 L=306)   936    158        3 ( 315/1263- 361/1416) -73.3 E57-E50 N05-N15
      3 1765 1334 (N29W54 L= 54)   444     80        0 (1731/1301-1797/1366) -26.6 W50-W58 N27-N31
      4  899 1081 (N02E08 L=351)  1043    400        2 ( 837/ 914- 982/1243)  81.9 E10-E05 S08-N13
      5 1639 1082 (N12W40 L= 39)   256     49        0 (1608/1064-1663/1099) -21.8 W37-W42 N10-N12
      6 1174  637 (S21W15 L= 14)  2141    613        1 ( 926/ 528-1423/ 741)  34.6 W00-W29 S31-S11
      7  617  594 (S32E23 L=336)   364     80        0 ( 575/ 581- 660/ 608)   8.1 E26-E19 S33-S30
  

    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