Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1554 (N29W22 L=355)   431    119        0 (1199/1513-1297/1580)  14.0 W18-W26 N26-N30
      2  812 1525 (N33E08 L=325)   221     33        0 ( 787/1509- 842/1538)   0.0 E09-E05 N31-N33
      3  370 1318 (N25E42 L=291)   257     43        0 ( 340/1304- 396/1334)  11.3 E44-E39 N23-N25
      4 1355 1294 (N11W23 L=356)   328     67        1 (1323/1273-1396/1320) -36.9 W21-W25 N09-N13
      5 1829 1292 (N06W61 L= 34)   525     44        0 (1806/1259-1846/1320) -76.0 W58-W62 N04-N08
      6  812 1280 (N17E10 L=323)   294     48        0 ( 783/1261- 845/1299)   0.0 E11-E07 N15-N17
      7 1855 1105 (S05W59 L= 32)  2877    321        3 (1746/ 982-1920/1221) -45.0 W50-W67 S12-N02
      8  213  982 (N06E57 L=276)  2188    284        6 ( 104/ 936- 333/1027)  -4.2 E72-E45 N02-N09
      9 1593  969 (S10W35 L=  8)   318     86        1 (1567/ 939-1630/1004) -26.6 W33-W37 S11-S07
     10 1895  837 (S21W65 L= 38)  3499    243        4 (1830/ 732-1930/ 948)  72.9 W57-W72 S27-S14
     11 1184  638 (S26W06 L=338)  3348    908        7 ( 758/ 517-1544/ 804)   7.5 E23-W30 S30-S19
     12  304  656 (S14E56 L=277)   499     92        2 ( 264/ 631- 351/ 678)  -7.1 E60-E51 S15-S12
     13 1810  658 (S32W59 L= 32)   387     33        2 (1788/ 643-1830/ 678) -15.9 W55-W62 S33-S30
     14  562  536 (S24E38 L=295)  1613    261        4 ( 491/ 474- 621/ 637)  34.7 E46-E32 S28-S18
     15 1631  516 (S40W42 L= 15)   533     69        1 (1598/ 492-1661/ 550)  53.1 W39-W45 S41-S37
     16 1398  440 (S42W21 L=354)   396     57        2 (1362/ 412-1425/ 464) -33.7 W17-W23 S44-S40
     17  584  343 (S38E47 L=286)  1187    203        1 ( 519/ 276- 639/ 427)  22.2 E60-E37 S41-S32
     18  792  221 (S51E38 L=295)  1155    131        1 ( 741/ 197- 858/ 255)   0.0 E45-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