Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  705 1785 (N39E51 L=140)  1425    163        1 ( 635/1735- 764/1839)  -8.4 E68-E41 N36-N41
      2  730 1640 (N29E37 L=154)   997    229        1 ( 617/1599- 815/1689)  50.7 E43-E32 N23-N34
      3 1567 1512 (N42W29 L=220)  1703    253        1 (1498/1407-1617/1602) -50.2 W22-W34 N35-N47
      4  413 1494 (N13E54 L=137)   633     88        0 ( 379/1454- 448/1533)  68.2 E55-E52 N10-N15
      5 1506 1341 (N28W24 L=215)  1843    499        1 (1355/1144-1583/1511) -43.5 W10-W32 N16-N35
      6  973 1321 (N14E10 L=181)   214     43        2 ( 952/1295- 993/1346) -33.7 E11-E08 N12-N15
      7 1808 1181 (N25W54 L=245)  1089    121        1 (1780/1119-1826/1237)  76.0 W50-W56 N21-N29
      8  277 1061 (S15E52 L=139)  2762    486        3 ( 170/ 888- 363/1219)  52.0 E66-E46 S27-S04
      9 1600  942 (N06W38 L=229)   280     69        0 (1565/ 920-1634/ 967)  45.0 W35-W39 N04-N08
     10  519  917 (S19E30 L=161)   273     85        0 ( 493/ 882- 540/ 957)  59.0 E31-E28 S21-S16
     11 1265  847 (S07W18 L=209)   212     72        0 (1228/ 829-1295/ 869)  45.0 W16-W19 S08-S05
     12 1490  800 (S05W34 L=225)   235     98        1 (1451/ 773-1537/ 818)  38.7 W31-W36 S07-S02
     13  892  717 (S23E01 L=190)   325    139        0 ( 836/ 680- 936/ 767)  66.8 E04-W00 S26-S19
     14 1439  730 (S10W32 L=223)   258     63        1 (1416/ 702-1458/ 762) -56.3 W30-W34 S11-S08
     15 1422  662 (S14W33 L=225)   451    105        2 (1384/ 622-1455/ 698) -26.6 W30-W36 S16-S12
     16  352  620 (S43E46 L=145)  1207    129        3 ( 284/ 590- 416/ 656)  11.9 E56-E37 S45-S39
     17 1318  629 (S19W28 L=219)   335     81        0 (1270/ 617-1354/ 652)  31.0 W24-W30 S20-S16
     18 1457  474 (S25W45 L=236)   672    237        1 (1368/ 401-1551/ 533)  59.7 W41-W49 S31-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