Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1664 1433 (N23W47 L=244)   641    114        1 (1620/1385-1711/1473) -29.7 W43-W50 N19-N25
      2  317 1327 (N17E51 L=146)   381     76        1 ( 300/1301- 345/1365) -33.7 E53-E47 N15-N19
      3  655 1311 (N15E23 L=173)   626    172        3 ( 579/1288- 739/1334)   0.0 E28-E17 N13-N16
      4  267 1232 (N11E54 L=143)   922    156        2 ( 249/1156- 289/1297) -83.7 E56-E51 N06-N15
      5  210 1226 (N11E60 L=136)   833    106        1 ( 184/1175- 235/1269)  59.0 E62-E57 N08-N13
      6  320 1169 (N07E48 L=149)   845    136        3 ( 283/1104- 357/1233)  54.5 E51-E45 N02-N10
      7  869  839 (S14E10 L=187)   843    328        0 ( 807/ 700- 902/ 989)  67.6 E14-E07 S22-S05
      8 1429  789 (S17W26 L=223)   256     90        0 (1390/ 774-1468/ 809)   0.0 W23-W28 S18-S16
      9  877  477 (S38E12 L=185)   807    206        0 ( 772/ 456- 974/ 511)   7.1 E20-E04 S39-S35
     10 1698  441 (S40W66 L=263)   916    114        0 (1648/ 424-1752/ 461)  -2.4 W56-W80 S40-S38
     11  911  322 (S50E11 L=186)   489    105        1 ( 864/ 303- 966/ 343)  10.3 E16-E05 S52-S48
  

    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