Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1817 1430 (N13W66 L= 68)  1227    256        1 (1763/1398-1878/1472) -19.7 W60-W75 N11-N16
      2  770 1280 (N18E12 L=350)   210     94        0 ( 749/1225- 788/1324) -90.0 E13-E11 N15-N21
      3  220 1264 (N25E55 L=307)   628    180        0 ( 174/1243- 266/1283) -21.8 E60-E50 N23-N27
      4  254 1135 (N17E51 L=311)  1466    194        3 ( 213/1058- 307/1226)  59.0 E55-E46 N12-N22
      5  136  979 (N09E67 L=296)  1908    284        7 (  96/ 914- 175/1028) -66.8 E73-E61 N05-N12
      6 1924  891 (S20W66 L= 68)  1730    252        6 (1884/ 831-1963/ 940)  56.3 W60-W72 S23-S16
      7  757  724 (S14E21 L=342)   435    189        0 ( 723/ 640- 791/ 818)  51.3 E24-E16 S18-S08
      8  618  732 (S12E29 L=333)   265     76        0 ( 603/ 690- 629/ 780)  63.4 E30-E27 S14-S08
      9  397  674 (S12E47 L=315)   374     59        0 ( 365/ 650- 428/ 696) -36.9 E49-E44 S13-S10
     10 1251  589 (S30W09 L= 11)   873    199        2 (1219/ 516-1281/ 695)  76.0 W05-W10 S34-S22
     11  285  539 (S19E63 L=299)  2210    395        4 ( 257/ 444- 333/ 649)  45.0 E69-E56 S24-S12
     12 1057  558 (S29E05 L=358)   242     70        0 (1024/ 527-1087/ 593)  21.8 E07-E02 S30-S26
  

    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