Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1334 (N30E19 L=  5)  1430    313        3 ( 529/1279- 837/1388)  -2.5 E30-E07 N27-N31
      2 1307 1218 (N13W20 L= 44)   629    147        2 (1244/1166-1379/1252) -36.9 W16-W24 N09-N15
      3  176 1043 (N16E64 L=320)   912    184        1 ( 139/1012- 220/1090) -24.0 E68-E58 N14-N18
      4 1667  985 (S07W41 L= 65)   863    232        1 (1622/ 872-1708/1085)  50.7 W36-W46 S12-S01
      5  674  724 (S07E26 L=357)  1640    307        4 ( 580/ 583- 720/ 881)  57.1 E35-E22 S15-N02
      6  936  793 (S06E09 L= 15)   644    196        2 ( 843/ 762-1036/ 824) -20.0 E14-E03 S09-S04
      7 1750  745 (S24W47 L= 71)  1118    143        2 (1699/ 702-1796/ 781)  -5.7 W42-W52 S26-S21
      8 1554  716 (S21W30 L= 54)   369    114        1 (1523/ 672-1574/ 769) -74.1 W28-W31 S24-S17
      9  879  514 (S24E18 L=  6)   251     52        1 ( 856/ 485- 905/ 539)  36.9 E19-E15 S25-S21
     10  648  305 (S37E43 L=340)  1335    218        2 ( 609/ 252- 685/ 376)  57.5 E48-E40 S41-S30
  

    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