Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  788 1703 (N30E33 L= 66)  1099    250        1 ( 679/1636- 869/1772) -15.3 E45-E23 N27-N34
      2  442 1654 (N22E58 L= 40)  3083    244        9 ( 359/1528- 522/1731) -56.3 E70-E49 N13-N28
      3 1168 1637 (N34E04 L= 94)   707    129        2 (1130/1585-1207/1678) -35.5 E07-W00 N30-N36
      4 1148 1260 (N09W02 L=101)   804    381        2 (1104/1095-1181/1423) -66.0 E02-W07 N00-N18
      5 1412 1347 (N20W17 L=115)   675    239        2 (1315/1273-1499/1415) -15.9 W09-W23 N16-N22
      6 1805 1275 (N26W50 L=149)   890    111        2 (1778/1221-1833/1336) -66.8 W47-W53 N23-N30
      7 1371  805 (S12W25 L=124)   237     61        1 (1350/ 774-1392/ 837) -36.9 W23-W27 S13-S10
      8 1477  766 (S12W33 L=131)   505    154        0 (1425/ 722-1551/ 801)  49.4 W30-W37 S15-S08
      9 1000  713 (S25W05 L=104)   346     86        1 ( 960/ 676-1043/ 746) -18.4 W02-W08 S26-S23
     10 1563  399 (S30W59 L=157)   978    214        2 (1472/ 352-1646/ 425)   6.8 W49-W74 S31-S27
     11  848  196 (S68W17 L=115)   810    120        2 ( 789/ 161- 896/ 229)  -5.0 W05-W28 S70-S65
  

    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