Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1290 1543 (N35W03 L= 77)  1425    227        2 (1232/1434-1336/1638) -47.5 E03-W08 N28-N40
      2  379 1490 (N09E57 L= 16)   794    138        0 ( 354/1447- 405/1550) -81.9 E59-E52 N06-N13
      3  677 1464 (N14E33 L= 41)   277     86        0 ( 656/1428- 710/1496)  78.7 E33-E31 N11-N16
      4 1599 1168 (N19W32 L=106)  1040    190        3 (1551/1081-1639/1256)  80.5 W29-W34 N13-N25
      5 1825 1039 (N18W55 L=129)  1398    163        1 (1794/ 957-1851/1117) -52.1 W50-W59 N13-N22
      6 1302 1040 (N04W15 L= 89)   347     75        0 (1258/1019-1336/1065)  45.0 W13-W17 N02-N06
      7 1902  868 (N11W71 L=144)   654     69        0 (1871/ 849-1922/ 900)  38.7 W65-W75 N09-N13
      8  955  743 (S21W03 L= 76)   494    230        0 ( 853/ 706-1068/ 786)  42.0 E02-W08 S25-S16
      9  539  712 (S33E25 L= 49)  1040    246        1 ( 448/ 677- 631/ 751)  11.3 E32-E17 S35-S29
     10 1595  568 (S14W51 L=124)   594    116        1 (1555/ 522-1636/ 608)  71.6 W48-W52 S17-S11
     11 1546  400 (S25W58 L=132)   570    112        2 (1493/ 363-1590/ 436)  81.9 W53-W62 S28-S21
  

    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