Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  954 1559 (N32E00 L=191)   811    190        1 ( 847/1530-1046/1588) -18.4 E07-W05 N29-N35
      2 1729 1356 (N14W52 L=243)   327     37        1 (1713/1326-1741/1385)  76.0 W50-W53 N11-N15
      3 1817 1180 (N03W57 L=248)   444     47        0 (1792/1158-1839/1204)  21.8 W54-W59 N01-N04
      4  120 1064 (N08E69 L=122)   607    179        1 (  99/1022- 140/1104) -59.0 E73-E65 N05-N10
      5  223 1001 (N03E56 L=135)   339     42        0 ( 205/ 977- 245/1022)  26.6 E58-E53 N02-N04
      6 1884  936 (S12W63 L=254)   521     64        0 (1861/ 903-1903/ 960) -76.0 W59-W65 S14-S10
      7  650  883 (S07E24 L=167)   668    146        1 ( 566/ 825- 714/ 931) -36.9 E29-E20 S11-S03
      8 1790  762 (S22W53 L=244)  3409    571        5 (1666/ 671-1926/ 884)  26.6 W42-W68 S27-S15
      9 1796  789 (S21W54 L=245)   377     83        1 (1779/ 753-1823/ 817)  76.0 W51-W56 S23-S19
     10  688  643 (S22E25 L=166)   225     93        0 ( 668/ 604- 719/ 683)  51.3 E26-E22 S24-S19
     11 1045  612 (S27E02 L=189)   392     50        0 (1013/ 590-1074/ 639)  45.0 E04-W00 S28-S25
     12 1637  409 (S46W49 L=240)   471    112        0 (1609/ 390-1666/ 432) -23.2 W45-W52 S48-S44
  

    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