Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  363 1203 (S10E47 L=231)   595    158        2 ( 340/1125- 381/1291) -73.3 E49-E44 S14-S04
      2 1653 1209 (N24W35 L=312)   283     83        1 (1639/1169-1666/1257) -78.7 W33-W35 N21-N26
      3  227 1138 (S16E57 L=220)   694    205        0 ( 174/1071- 268/1214)  68.2 E61-E53 S21-S11
      4  530  905 (S23E27 L=250)   982    263        3 ( 412/ 856- 635/ 955)   0.0 E36-E19 S24-S21
      5  192  872 (S33E57 L=221)  1557    396        1 ( 105/ 833- 264/ 912)   3.2 E67-E49 S35-S31
      6  978  840 (S16W02 L=280)  1073    382        4 ( 787/ 792-1138/ 893)  15.9 E10-W11 S18-S11
      7 1227  405 (S34W34 L=312)   257     71        0 (1189/ 394-1264/ 418)  36.9 W31-W35 S35-S32
  

    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