Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1095 1366 (N26W07 L=194)   407    128        1 (1043/1334-1148/1395) -29.7 W03-W10 N23-N27
      2  559 1341 (N27E31 L=156)   300     61        1 ( 524/1325- 590/1357)  -9.5 E34-E28 N26-N28
      3  356 1070 (N11E46 L=141)  1115    251        2 ( 331/ 956- 386/1191) -71.6 E48-E43 N03-N18
      4 1652  941 (S06W41 L=229)   917    156        2 (1585/ 899-1729/ 969)   0.0 W36-W48 S07-S03
      5 1851  891 (S11W62 L=249)   561     97        1 (1821/ 860-1880/ 935) -35.5 W57-W66 S13-S08
      6 1551  760 (S16W33 L=221)   235     92        0 (1507/ 729-1578/ 784) -31.0 W30-W35 S18-S13
      7  486  746 (S10E38 L=149)   358     51        1 ( 462/ 722- 515/ 770) -45.0 E39-E35 S11-S08
      8  969  546 (S25E07 L=180)  2048    561        3 ( 830/ 342-1148/ 710)  43.9 E21-W05 S40-S15
      9 1650  618 (S27W44 L=231)   684    178        0 (1622/ 539-1685/ 700)  84.8 W42-W46 S32-S21
     10 1756  602 (S30W57 L=245)   573     60        1 (1720/ 582-1786/ 626)   7.1 W53-W61 S30-S28
     11 1468  577 (S28W28 L=216)   500    117        0 (1429/ 529-1518/ 619)  39.8 W25-W31 S30-S24
  

    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