Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  584 1513 (N48E18 L= 41)   529    103        1 ( 541/1486- 635/1529) -12.5 E22-E13 N46-N50
      2  889 1102 (N15E06 L= 54)   184     52        1 ( 861/1083- 917/1119) -14.0 E07-E03 N13-N15
      3  259  873 (N15E56 L=  3)  3127    450        4 ( 187/ 728- 393/1057) -61.4 E62-E47 N04-N27
      4 1767  716 (S34W44 L=103)   582     88        2 (1732/ 682-1798/ 750)   8.1 W39-W47 S35-S31
      5 1257  540 (S28W00 L= 60)   293     72        1 (1220/ 515-1292/ 560)   0.0 E02-W03 S28-S26
      6  455  413 (S16E58 L=  1)  1341    203        1 ( 385/ 352- 533/ 489) -59.0 E63-E53 S20-S09
      7 1089  360 (S35E17 L= 43)   989    209        2 (1024/ 298-1144/ 448)  25.0 E24-E09 S38-S30
      8  888  179 (S44E45 L= 15)  1341    170        2 ( 853/ 111- 936/ 241)  20.1 E64-E34 S50-S39
  

    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