Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  558 1662 (N18E50 L=225)   406     38        0 ( 533/1642- 584/1684)  63.4 E52-E48 N16-N20
      2  345 1423 (S00E56 L=219)   493     56        0 ( 329/1390- 361/1454) -36.9 E58-E53 S02-N01
      3  808 1381 (N08E22 L=253)   583    164        1 ( 749/1330- 881/1448) -15.3 E27-E16 N05-N10
      4 1275 1222 (N11W08 L=284)  1334    325        0 (1238/1043-1301/1370) -60.9 W02-W13 N01-N19
      5  187 1193 (S16E63 L=212)  1043    134        1 ( 157/1142- 221/1235)  80.5 E66-E59 S19-S13
      6 1485  958 (N03W28 L=303)   387    109        0 (1454/ 910-1504/1021) -71.6 W25-W29 N00-N06
      7  355  944 (S28E41 L=234)  1414    384        1 ( 220/ 893- 484/1019)  17.7 E54-E31 S32-S20
      8  968  896 (S15W00 L=276)  1230    367        0 ( 774/ 854-1136/ 941)  15.9 E11-W10 S18-S11
  

    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