Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1465 1650 (N29W52 L=251)   506     80        0 (1424/1625-1502/1674) -68.2 W49-W54 N26-N31
      2  683 1581 (N50E08 L=191)   253     44        1 ( 660/1562- 712/1605) -11.3 E10-E05 N48-N51
      3 1564 1536 (N20W51 L=250)   378     67        0 (1527/1524-1606/1550) -15.9 W48-W55 N18-N20
      4 1533 1391 (N13W41 L=240)   341     62        1 (1512/1360-1558/1420)  21.8 W39-W44 N11-N14
      5  198 1378 (N47E69 L=129)   795     88        0 ( 162/1347- 241/1397)  -7.1 E82-E58 N44-N48
      6 1596 1355 (N09W44 L=243)   313     70        1 (1563/1329-1622/1386) -11.3 W41-W46 N08-N11
      7  383 1130 (N28E40 L=159)   252     42        0 ( 354/1116- 407/1147) -11.3 E42-E37 N26-N28
      8 1040 1079 (N10W02 L=201)   369    114        1 ( 996/1030-1084/1118) -71.6 W00-W03 N06-N12
      9  762 1082 (N16E14 L=185)   288    140        0 ( 707/1048- 808/1108) -45.0 E17-E11 N13-N18
     10 1632  903 (S16W34 L=232)   262     49        1 (1610/ 881-1647/ 935)  45.0 W31-W35 S17-S14
     11 1199  679 (S17W01 L=200)   898    256        3 (1083/ 617-1280/ 749)   4.4 E06-W07 S21-S14
     12  340  609 (S02E57 L=142)   514     87        0 ( 298/ 579- 385/ 632) -33.7 E60-E53 S04-N00
     13  252  572 (S04E70 L=129)   543     76        2 ( 226/ 538- 282/ 605) -53.1 E71-E67 S05-S01
     14 1616  455 (S47W28 L=227)  1433    185        2 (1543/ 418-1695/ 525)   3.6 W20-W37 S49-S43
     15  385  445 (S13E62 L=137)  1143    132        2 ( 364/ 394- 410/ 507)  45.0 E68-E58 S16-S08
     16  751  332 (S28E39 L=160)  2537    431        2 ( 714/ 149- 799/ 502)  39.1 E59-E27 S43-S17
     17  478  425 (S16E54 L=144)   945    148        5 ( 445/ 382- 536/ 481) -76.0 E57-E50 S19-S11
     18  487  312 (S24E64 L=135)   990    107        3 ( 431/ 290- 534/ 332) -18.4 E70-E58 S26-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