Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1749 (N44E08 L= 62)   591     79        1 ( 965/1716-1045/1768)  31.0 E10-E04 N40-N45
      2 1374 1601 (N35W21 L= 91)  4004    734        8 (1116/1463-1674/1714) -14.7 W02-W45 N27-N41
      3  752 1624 (N32E23 L= 46)   741    127        0 ( 708/1567- 791/1673) -41.2 E27-E19 N28-N35
      4 1261 1280 (N12W13 L= 83)   584    168        1 (1238/1185-1292/1364) -74.7 W11-W15 N06-N17
      5  614 1310 (N11E27 L= 43)   343    106        0 ( 563/1282- 652/1329) -18.4 E30-E24 N08-N11
      6  938 1243 (N08E06 L= 63)   257     66        0 ( 913/1211- 964/1283) -53.1 E08-E04 N06-N10
      7 1833 1218 (N14W56 L=126)   588    125        0 (1797/1193-1880/1251) -14.0 W52-W61 N12-N15
      8  306 1202 (N03E49 L= 20)   378     88        0 ( 272/1163- 336/1239)  38.7 E52-E47 N01-N05
      9  689  707 (S26E20 L= 50)   295     80        0 ( 655/ 672- 722/ 745)  45.0 E22-E17 S27-S22
     10  271  582 (S34E57 L= 13)   725     92        0 ( 247/ 553- 309/ 622) -68.2 E59-E54 S36-S31
     11 1645  511 (S31W51 L=121)   864    283        0 (1546/ 471-1766/ 566)  25.2 W43-W60 S34-S25
  

    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