Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1361 (N33E27 L= 32)   327     99        0 ( 532/1344- 627/1382) -15.9 E30-E23 N31-N34
      2 1636 1313 (N13W46 L=105)   350    145        1 (1591/1276-1667/1358)  18.4 W41-W50 N11-N15
      3  298 1279 (N29E51 L=  8)  2014    345        6 ( 201/1226- 418/1345)   5.7 E60-E40 N26-N33
      4  752 1101 (N14E16 L= 43)   362    105        0 ( 704/1076- 806/1121) -26.6 E19-E12 N12-N15
      5  247  664 (S09E64 L=355)  2796    424        8 ( 209/ 545- 279/ 804)  67.6 E74-E59 S17-N00
      6  385  701 (S07E48 L= 11)  1007    223        2 ( 343/ 636- 461/ 741)  15.3 E53-E41 S10-S04
      7 1269  615 (S23W10 L= 69)   222     90        0 (1225/ 586-1311/ 649)  26.6 W06-W13 S24-S20
      8  368  590 (S14E53 L=  6)   339     52        0 ( 351/ 566- 384/ 616) -90.0 E55-E51 S15-S12
  

    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