Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  478 1622 (N20E51 L=220)   519    154        1 ( 434/1602- 530/1649)   0.0 E55-E46 N18-N21
      2  324 1232 (S06E50 L=221)   385    107        1 ( 307/1188- 339/1276) -68.2 E51-E47 S08-S03
      3  210 1175 (S11E60 L=211)  2562    316        8 ( 166/1066- 248/1275) -55.3 E66-E55 S17-S04
      4  370 1088 (S14E43 L=228)   588    219        5 ( 302/1063- 457/1119)  14.0 E48-E36 S15-S12
      5  163 1088 (S16E63 L=208)   512    150        4 ( 132/1065- 195/1112)  15.9 E67-E59 S17-S14
      6  932 1035 (S08E06 L=266)   284     53        1 ( 907/1006- 956/1072)  78.7 E06-E04 S10-S05
      7  406 1037 (S17E40 L=232)   262     45        0 ( 386/1011- 428/1065)  63.4 E41-E38 S18-S14
      8  791  979 (S14E13 L=258)   359    129        1 ( 710/ 960- 847/ 994)  20.6 E17-E09 S15-S12
      9 1294  846 (S11W19 L=290)   387    129        2 (1260/ 812-1345/ 892) -18.4 W17-W23 S13-S08
     10 1219  790 (S16W16 L=288)   746    284        0 (1101/ 749-1366/ 867)   0.0 W07-W25 S18-S14
     11  830  761 (S26E07 L=265)  1456    472        4 ( 639/ 659-1053/ 849)  39.6 E18-W05 S35-S16
     12 1211  593 (S28W22 L=293)   251     73        0 (1175/ 580-1252/ 615)  11.3 W19-W24 S29-S26
     13 1836  578 (S11W74 L=345)   724    146        2 (1804/ 546-1858/ 603)  53.1 W68-W77 S13-S09
     14  620  539 (S45E19 L=252)   571    128        1 ( 600/ 486- 651/ 596) -82.9 E20-E16 S48-S40
     15 1579  554 (S21W48 L=319)   515     62        0 (1546/ 523-1618/ 586)  59.0 W45-W50 S23-S18
     16  951  498 (S40W08 L=279)   428    146        1 ( 881/ 466-1031/ 533)  41.6 W03-W12 S43-S35
     17  741  448 (S48E07 L=264)   270     46        1 ( 715/ 431- 764/ 473)  14.0 E09-E05 S49-S46
     18 1211  397 (S40W30 L=301)   384     55        1 (1177/ 377-1247/ 419)  45.0 W27-W32 S41-S37
  

    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