Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1314 1764 (N53W40 L= 54)   835    109        1 (1276/1728-1364/1783)   0.0 W33-W47 N50-N55
      2  553 1462 (N35E34 L=340)   713    210        0 ( 516/1374- 601/1534)  70.0 E36-E31 N28-N40
      3  717 1432 (N32E20 L=355)   916    229        2 ( 648/1374- 798/1518) -34.7 E26-E13 N28-N38
      4  739 1218 (N18E18 L=357)   375    145        0 ( 688/1164- 797/1260)  35.5 E20-E13 N14-N20
      5 1262 1159 (N12W16 L= 30)   209     35        0 (1232/1144-1285/1174)   0.0 W13-W17 N10-N12
      6 1262 1090 (N07W15 L= 30)   312     60        1 (1228/1064-1293/1109) -26.6 W13-W17 N05-N08
      7 1190 1021 (N04W10 L= 25)   188     41        0 (1169/ 998-1209/1041)  45.0 W08-W11 N02-N04
      8 1605  816 (S12W38 L= 52)   478    123        1 (1571/ 766-1653/ 845)  33.7 W35-W41 S15-S10
      9  329  520 (S26E62 L=313)  1166    190        0 ( 272/ 436- 373/ 588) -90.0 E65-E58 S32-S21
  

    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