Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1582 (N11E62 L=206)   680     78        1 ( 333/1552- 397/1608)  36.9 E66-E58 N09-N13
      2 1696 1414 (N35W35 L=304)  1495    265        0 (1575/1337-1790/1473)  -2.7 W24-W45 N32-N36
      3  324 1346 (S04E53 L=215)   524    115        1 ( 292/1304- 379/1390)  50.2 E56-E49 S06-N00
      4  333 1231 (S10E49 L=219)   481     88        1 ( 305/1208- 371/1263) -15.9 E52-E45 S11-S08
      5  943 1189 (S00E08 L=260)   188     47        0 ( 929/1160- 956/1226) -45.0 E09-E07 S01-N01
      6 1858 1061 (N18W55 L=323)   441     39        1 (1837/1035-1880/1090)  76.0 W52-W58 N16-N20
      7 1194  661 (S23W20 L=288)   633    114        0 (1153/ 623-1244/ 710)  -7.1 W15-W23 S24-S21
      8 1167  605 (S27W20 L=288)   231     56        0 (1143/ 577-1187/ 638) -21.8 W17-W22 S28-S25
      9  713  495 (S46E06 L=262)  1622    394        0 ( 560/ 432- 895/ 568)   3.9 E20-W09 S48-S43
  

    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