Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  490 1565 (N53E36 L= 73)   517     51        0 ( 472/1527- 513/1595) -56.3 E39-E32 N49-N55
      2  730 1566 (N47E10 L= 99)   766    120        2 ( 680/1534- 786/1594) -36.9 E14-E04 N43-N50
      3 1723 1465 (N15W63 L=171)   455     86        0 (1694/1428-1742/1498) -90.0 W60-W65 N12-N17
      4 1799 1452 (N12W72 L=181)   994     86        3 (1779/1410-1817/1491)  59.0 W68-W79 N09-N14
      5  704 1294 (N29E15 L= 93)   360    134        0 ( 666/1245- 756/1338)  32.0 E18-E10 N26-N31
      6  192 1112 (N25E64 L= 45)  1961    321        4 ( 124/1000- 249/1209)  33.7 E75-E56 N18-N31
      7  336 1118 (N24E47 L= 62)   765    155        1 ( 290/1058- 374/1182) -45.0 E51-E43 N20-N28
      8  807  828 (S00E17 L= 92)   657    192        2 ( 755/ 740- 860/ 893) -73.3 E18-E14 S06-N04
      9 1083  749 (S10E02 L=107)   248     58        1 (1055/ 717-1113/ 776)  14.0 E04-W00 S11-S08
     10 1723  721 (S28W43 L=152)   482     49        1 (1700/ 692-1743/ 752)  71.6 W40-W45 S29-S26
     11  468  385 (S24E59 L= 50)  1888    274        1 ( 382/ 308- 550/ 438) -57.5 E65-E53 S30-S19
  

    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