Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  735 1209 (N24E12 L= 86)  1586    479        2 ( 528/1091- 912/1324) -47.7 E24-E03 N13-N36
      2 1931 1192 (S13W72 L=171)  3313    167        7 (1898/1103-1965/1267)  84.3 W66-W81 S19-S09
      3 1747 1050 (S14W45 L=144)  1458    443        2 (1620/ 923-1784/1230)  31.0 W32-W52 S20-S04
      4 1866  961 (S23W56 L=154)  1131    211        1 (1805/ 871-1920/1024)  16.4 W47-W64 S26-S20
      5 1337  741 (S18W10 L=108)   929    314        0 (1283/ 602-1395/ 871)  42.7 W03-W16 S24-S12
      6 1524  792 (S21W23 L=121)   263    103        0 (1503/ 738-1536/ 840)  51.3 W20-W24 S23-S18
      7  911  786 (S04E12 L= 86)   188     63        1 ( 889/ 757- 932/ 813) -76.0 E13-E11 S05-S01
      8 1053  715 (S11E07 L= 92)   282    104        1 (1036/ 657-1075/ 757)  45.0 E09-E04 S14-S09
      9  597  527 (S11E40 L= 58)   698    204        3 ( 567/ 426- 619/ 612)  57.5 E45-E36 S16-S05
     10  291  535 (S05E67 L= 32)  1321    109        2 ( 256/ 490- 332/ 572)  45.0 E72-E61 S07-S02
     11 1348  545 (S30W06 L=104)   283     95        0 (1294/ 531-1398/ 559) -33.7 W02-W08 S32-S28
     12  462  387 (S17E59 L= 39)  1422    166        1 ( 424/ 317- 513/ 465)  70.0 E62-E55 S22-S11
  

    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