Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1209 1698 (N38W35 L=287)  1255    230        2 (1105/1666-1333/1718) -22.4 W27-W44 N34-N41
      2 1687 1649 (N21W73 L=325)  2202    124        1 (1638/1609-1746/1685) -37.9 W67-W81 N16-N23
      3 1467 1486 (N18W41 L=292)  1814    391        1 (1356/1366-1607/1619) -76.0 W36-W46 N08-N28
      4 1602 1523 (N17W53 L=304)  2353    467        2 (1435/1392-1766/1614) -52.7 W45-W61 N05-N26
      5  605 1306 (N30E20 L=232)   745    245        0 ( 510/1277- 713/1334) -29.7 E26-E12 N25-N34
      6  288 1306 (N38E49 L=202)  1899    308        1 ( 137/1286- 433/1326)  -6.3 E78-E33 N34-N41
      7  922 1118 (N12E04 L=248)   305     89        0 ( 894/1074- 947/1159) -90.0 E04-E02 N08-N14
      8 1684 1093 (S09W42 L=294)   515    123        1 (1631/1043-1733/1135)   5.7 W37-W47 S10-S07
      9 1914  895 (S28W66 L=317)  1159     81        1 (1869/ 874-1956/ 917)  -9.0 W58-W77 S29-S25
     10 1752  807 (S28W44 L=295)   761    187        1 (1706/ 723-1795/ 882)  37.9 W39-W48 S31-S23
     11 1618  663 (S33W30 L=282)   668    144        0 (1578/ 622-1667/ 739) -61.4 W27-W35 S37-S26
     12 1490  660 (S29W20 L=271)   358    136        1 (1444/ 611-1537/ 719)  32.0 W16-W24 S31-S26
     13  580  644 (S08E36 L=215)   564    177        0 ( 528/ 568- 607/ 714) -90.0 E37-E34 S12-S02
     14  495  595 (S09E44 L=207)   324     46        0 ( 472/ 572- 518/ 619) -90.0 E45-E42 S10-S07
     15  847  518 (S21E24 L=227)   635    131        1 ( 822/ 459- 875/ 591)  48.8 E27-E19 S24-S16
     16 1324  412 (S42W03 L=254)   680    152        0 (1258/ 365-1374/ 461)  16.7 E03-W07 S43-S39
     17 1146  285 (S46E16 L=236)   312     59        0 (1116/ 269-1182/ 304)   8.1 E18-E11 S47-S44
     18  951  167 (S50E41 L=210)  1460    172        3 ( 873/ 128-1026/ 205)   4.6 E55-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