Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  503 1672 (N20E56 L=314)  1900    246        3 ( 445/1628- 574/1737)  66.0 E61-E53 N16-N25
      2 1334 1591 (N38W03 L= 13)   511     81        1 (1297/1556-1361/1618)   0.0 W00-W04 N35-N40
      3  761 1518 (N16E30 L=340)   533    163        0 ( 738/1460- 805/1589) -36.9 E33-E25 N13-N20
      4 1427 1463 (N32W12 L= 23)  1649    269        4 (1371/1354-1472/1567) -58.4 W08-W17 N26-N39
      5  492 1465 (N07E46 L=325)   515    112        2 ( 462/1425- 530/1501)  59.0 E48-E43 N03-N08
      6  587 1424 (N06E37 L=333)   236     92        0 ( 550/1388- 606/1465)  80.5 E38-E35 N03-N09
      7 1799 1186 (N28W47 L= 57)   494     87        0 (1773/1134-1815/1229) -51.3 W44-W49 N24-N29
      8 1393 1068 (N08W20 L= 30)   197     45        0 (1376/1044-1410/1097) -76.0 W18-W20 N05-N09
      9 1743 1027 (N16W45 L= 55)   416     85        2 (1718/1002-1780/1051)   9.5 W42-W48 N14-N17
     10 1905  926 (N16W65 L= 75)  2259    247        6 (1881/ 823-1927/1019) -54.0 W60-W72 N11-N22
     11  897  925 (S14E04 L=  6)   248    146        1 ( 823/ 896- 966/ 968)  48.8 E08-E01 S17-S09
     12  255  892 (S32E50 L=320)   854    216        1 ( 224/ 807- 315/ 944)  55.0 E53-E44 S37-S27
     13  303  765 (S39E45 L=326)   577    109        1 ( 271/ 736- 335/ 804) -18.4 E47-E41 S40-S36
     14 1431  698 (S11W33 L= 43)   604    174        0 (1359/ 634-1514/ 736)  56.3 W30-W37 S16-S07
     15 1144  708 (S19W16 L= 26)   380    153        0 (1066/ 686-1224/ 729)  16.7 W11-W21 S21-S16
     16 1197  574 (S25W24 L= 34)   610    169        0 (1114/ 562-1281/ 595)  33.7 W19-W28 S27-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