Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1700 (N38W36 L=287)  1236    246        1 (1108/1668-1335/1718) -23.6 W28-W44 N34-N41
      2 1687 1653 (N21W74 L=325)  2130    133        1 (1641/1611-1744/1691) -58.0 W67-W81 N16-N24
      3 1472 1489 (N18W41 L=292)  1647    393        0 (1362/1368-1609/1623) -78.7 W37-W46 N08-N28
      4 1575 1548 (N19W52 L=303)  1853    339        0 (1437/1481-1693/1618) -51.3 W45-W59 N11-N26
      5  607 1308 (N30E19 L=231)   669    258        0 ( 514/1276- 719/1334) -29.7 E26-E12 N25-N33
      6  295 1306 (N38E48 L=203)  1748    312        0 ( 136/1283- 425/1327)  -7.4 E80-E34 N34-N41
      7  925 1127 (N12E03 L=248)   286     89        0 ( 897/1086- 954/1171) -78.7 E04-E02 N09-N14
      8 1687 1094 (S09W43 L=294)   587    136        2 (1636/1033-1739/1141)  24.0 W37-W47 S11-S07
      9 1913  898 (S28W66 L=317)  1008     88        1 (1872/ 875-1960/ 916)  -5.2 W58-W80 S28-S25
     10 1755  808 (S28W44 L=295)   642    174        0 (1711/ 732-1795/ 882)  41.2 W39-W49 S31-S23
     11 1624  667 (S33W31 L=282)   630    143        0 (1579/ 625-1675/ 740) -61.4 W28-W36 S37-S26
     12  586  644 (S08E36 L=215)   558    172        0 ( 533/ 568- 612/ 718)  84.3 E37-E34 S12-S02
     13  501  597 (S09E44 L=207)   267     44        1 ( 481/ 576- 520/ 620) -90.0 E45-E42 S10-S07
     14  851  517 (S21E24 L=227)   606    126        1 ( 821/ 460- 877/ 592)  41.2 E27-E19 S24-S16
     15 1307  402 (S42W01 L=252)   850    208        1 (1215/ 339-1373/ 475)  19.7 E07-W07 S44-S38
     16  785  205 (S42E49 L=201)   510     56        1 ( 755/ 183- 812/ 237)  16.7 E54-E44 S43-S39
     17  950  168 (S50E41 L=209)  1464    180        4 ( 861/ 130-1027/ 205)   0.0 E54-E30 S52-S47
  

    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