Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  561 1129 (N23E27 L=318)   579    193        0 ( 477/1102- 651/1157) -26.6 E32-E20 N20-N26
      2 1560  994 (S08W31 L= 16)   460    127        0 (1522/ 948-1618/1018)   7.1 W27-W35 S09-S06
      3  869  869 (N02E12 L=333)   727    194        1 ( 849/ 770- 894/ 961)  68.2 E15-E10 S03-N07
      4 1703  827 (S22W39 L= 24)   834    234        0 (1649/ 712-1723/ 942)  58.4 W33-W41 S28-S15
      5 1405  847 (S12W18 L=  3)   739    185        1 (1346/ 760-1437/ 932)  52.1 W12-W21 S16-S07
      6 1257  718 (S16W06 L=351)   325     95        0 (1212/ 694-1309/ 740)   0.0 W03-W09 S17-S14
      7 1371  594 (S27W11 L=356)  1133    326        1 (1210/ 473-1497/ 719)  19.2 E02-W21 S30-S22
      8 1728  666 (S34W40 L= 25)   565     75        0 (1690/ 635-1765/ 700)   8.1 W36-W43 S35-S32
      9 1027  493 (S24E13 L=332)   310     64        1 ( 999/ 466-1065/ 524)  21.8 E15-E10 S25-S23
     10  968  272 (S39E27 L=318)   473    122        1 ( 946/ 231- 995/ 315)  33.7 E31-E22 S42-S36
  

    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