Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  948 1542 (N26E09 L=197)   775    247        0 ( 847/1467-1036/1593) -21.0 E15-E02 N22-N29
      2  287 1370 (N13E55 L=151)  3001    336        8 ( 213/1236- 381/1522) -43.5 E67-E43 N04-N22
      3  509 1436 (N17E37 L=169)  1201    228        4 ( 428/1376- 585/1519)  45.0 E43-E32 N13-N22
      4  374 1406 (N15E47 L=159)   506     71        1 ( 350/1376- 402/1444)  90.0 E49-E44 N13-N17
      5 1427 1282 (N14W23 L=229)   777    276        0 (1287/1255-1527/1308)  14.9 W14-W29 N10-N15
      6  203 1016 (S09E58 L=148)  1357    309        7 ( 167/ 899- 263/1115)  77.0 E63-E52 S15-S02
      7  622 1041 (S06E25 L=182)   263     89        1 ( 599/1003- 641/1079) -63.4 E26-E23 S08-S04
      8  536  974 (S10E30 L=176)   263     59        0 ( 516/ 947- 562/ 996)  71.6 E31-E28 S12-S09
      9  630  933 (S13E24 L=182)   386     90        0 ( 591/ 898- 667/ 970)  45.0 E26-E21 S14-S10
     10  390  891 (S16E41 L=165)  1569    316        3 ( 318/ 804- 468/ 948) -41.2 E47-E35 S21-S12
     11 1177  927 (S09W10 L=216)   193     66        0 (1142/ 903-1202/ 941) -14.0 W07-W11 S10-S08
     12  205  856 (S18E60 L=147)  1850    169        6 ( 149/ 792- 280/ 909)  23.2 E67-E51 S21-S14
     13  538  835 (S19E30 L=176)   371    108        1 ( 500/ 803- 594/ 870) -26.6 E32-E26 S20-S16
     14  217  674 (S29E64 L=142)   695     60        1 ( 197/ 644- 236/ 722)  78.7 E67-E61 S31-S26
     15 1348  504 (S34W28 L=234)   566    154        2 (1273/ 462-1407/ 536)  -9.5 W21-W33 S36-S32
     16  792  506 (S39E13 L=193)   301    107        0 ( 740/ 490- 848/ 521)  -6.3 E17-E08 S39-S37
     17  965  450 (S41W00 L=207)   370    135        2 ( 905/ 434-1041/ 473)  16.7 E04-W06 S42-S39
     18 1416  240 (S52W54 L=260)  1381    189        4 (1336/ 222-1505/ 265)   2.1 W42-W70 S54-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