Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1286 1803 (N51W51 L=171)  1567    225        2 (1179/1759-1381/1843)  -1.9 W35-W66 N48-N52
      2 1525 1690 (N36W61 L=181)   431     49        1 (1503/1669-1554/1706) -26.6 W57-W64 N34-N37
      3 1728 1326 (N10W55 L=175)  1302    333        0 (1596/1275-1865/1405)  -3.0 W42-W80 N08-N11
      4  899 1227 (N21E04 L=116)   204    111        0 ( 876/1177- 935/1267)  38.7 E06-E01 N17-N22
      5 1048 1220 (N18W05 L=125)   256    102        0 ( 992/1197-1090/1236) -26.6 W01-W07 N15-N19
      6  637 1072 (N15E23 L= 96)   583    182        3 ( 579/1008- 697/1145) -63.4 E26-E20 N10-N20
      7 1780  932 (S14W51 L=170)  3284    588        9 (1628/ 654-1864/1100)  40.2 W35-W62 S28-S05
      8 1709 1035 (S06W45 L=165)   343     51        1 (1690/1005-1727/1062)  71.6 W43-W46 S07-S04
      9 1730  973 (S10W46 L=166)   603     81        1 (1702/ 945-1761/1014)  45.0 W43-W48 S12-S08
     10  697  848 (N01E23 L= 97)   246     70        1 ( 680/ 811- 717/ 879)  56.3 E23-E21 S01-N02
     11 1022  823 (S05E04 L=116)   189     86        0 ( 999/ 790-1044/ 861) -78.7 E04-E02 S07-S02
     12 1157  746 (S12W03 L=123)   531    124        5 (1109/ 700-1186/ 808)  56.3 W00-W05 S14-S08
     13  368  534 (S15E57 L= 62)  1817    381        1 ( 307/ 399- 453/ 698)  84.8 E63-E51 S25-S03
     14  851  601 (S16E18 L=102)   455     96        3 ( 804/ 572- 897/ 646)  33.7 E21-E15 S18-S13
     15  911  516 (S23E17 L=103)   425    104        1 ( 890/ 459- 941/ 555)  56.3 E19-E14 S26-S20
     16 1037  428 (S31E11 L=109)   695    121        2 ( 980/ 401-1097/ 458)   0.0 E15-E06 S32-S29
     17  623  403 (S27E41 L= 79)   862     75        3 ( 594/ 355- 659/ 442)  60.3 E44-E37 S30-S23
     18 1287  268 (S52W05 L=125)   677    100        1 (1245/ 245-1321/ 297)  20.6 W00-W08 S54-S49
  

    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