Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1784 (N31E65 L= 75)  1673    134        1 ( 543/1715- 597/1846) -29.7 E79-E57 N25-N37
      2 1508 1695 (N56W19 L=159)  1428    118        4 (1457/1658-1551/1744)  47.7 W13-W24 N49-N61
      3 1219 1654 (N40E06 L=134)   488    121        1 (1158/1620-1263/1671)   0.0 E10-E01 N38-N41
      4  696 1563 (N19E36 L=105)   283     19        0 ( 679/1539- 714/1586) -33.7 E37-E34 N17-N21
      5 1622 1570 (N49W32 L=172)   410     33        0 (1600/1550-1648/1586)  18.4 W29-W35 N46-N51
      6  918 1461 (N19E18 L=123)  1318    252        2 ( 864/1364- 975/1583) -38.2 E25-E11 N13-N25
      7  394 1516 (N10E57 L= 83)  2331    286        6 ( 315/1456- 451/1581)  41.2 E65-E50 N05-N14
      8 1346 1522 (N35W07 L=147)   343     54        0 (1313/1498-1377/1539)   0.0 W04-W09 N33-N36
      9 1063 1269 (N11E04 L=137)   206     23        0 (1048/1246-1084/1288) -45.0 E05-E02 N09-N12
     10 1674 1160 (N21W38 L=178)   465     90        0 (1623/1140-1717/1176)   6.3 W33-W42 N20-N22
     11  400  834 (S30E37 L=103)   362     95        0 ( 371/ 811- 442/ 870)  -8.1 E40-E33 S31-S28
     12  727  803 (S24E12 L=128)   381    129        0 ( 676/ 786- 787/ 833)  14.0 E16-E08 S25-S21
     13  802  758 (S25E06 L=134)   265    117        0 ( 745/ 737- 849/ 785)  45.0 E09-E04 S26-S21
     14  667  745 (S29E15 L=125)   309     48        1 ( 647/ 715- 695/ 775)  90.0 E16-E13 S31-S26
     15 1652  584 (S11W55 L=195)   414     44        0 (1634/ 555-1671/ 609) -76.0 W52-W56 S13-S09
     16 1434  567 (S19W39 L=179)   437     68        0 (1411/ 533-1455/ 600) -36.9 W36-W41 S19-S16
     17  889  526 (S37W06 L=146)   406     44        0 ( 860/ 508- 926/ 550)  21.8 W03-W08 S38-S35
     18 1494  421 (S25W52 L=193)   518     75        0 (1457/ 404-1539/ 439)   6.3 W48-W57 S26-S24
     19  968  281 (S52W26 L=166)   552     63        0 ( 941/ 253- 998/ 315) -18.4 W21-W30 S53-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