Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  814 1789 (N60W16 L=191)   351     64        0 ( 786/1765- 843/1811)   7.1 W11-W20 N58-N62
      2  512 1245 (N30E27 L=148)   381     68        1 ( 483/1213- 544/1281) -68.2 E28-E24 N28-N33
      3  510 1133 (N24E29 L=146)   371     82        1 ( 494/1089- 536/1172)  53.1 E30-E26 N21-N25
      4  191  981 (N21E61 L=114)   938    131        2 ( 171/ 925- 209/1045)  66.8 E64-E58 N18-N25
      5 1407  757 (S20W15 L=190)   555    228        0 (1334/ 680-1497/ 825)  17.1 W09-W22 S23-S17
      6 1656  777 (S27W34 L=209)   337     92        0 (1631/ 735-1671/ 821)  45.0 W31-W35 S28-S24
      7 1716  644 (S38W39 L=214)   441     49        1 (1690/ 623-1745/ 667) -39.8 W36-W42 S40-S35
      8 1332  536 (S32W05 L=180)   443    157        1 (1256/ 505-1391/ 571)   0.0 W00-W10 S33-S30
      9  451  471 (S13E54 L=121)   426    101        0 ( 432/ 421- 470/ 515)  50.2 E56-E51 S16-S10
     10 1122  353 (S38E15 L=160)   425     99        0 (1074/ 332-1171/ 377)  -7.1 E18-E10 S39-S37
  

    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