Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  418 1588 (N16E55 L=177)   458    109        0 ( 382/1557- 466/1610)  38.7 E57-E51 N13-N17
      2  524 1575 (N16E45 L=186)   293     61        0 ( 496/1548- 553/1607)  63.4 E46-E44 N14-N18
      3  692 1086 (S10E21 L=211)   952    315        3 ( 608/ 986- 743/1158) -35.0 E26-E16 S15-S05
      4  507 1049 (S15E32 L=199)   232     92        1 ( 469/1021- 549/1072)  -8.1 E35-E28 S16-S13
      5  623  885 (S23E22 L=210)   458    184        0 ( 538/ 851- 687/ 946)  37.9 E27-E18 S25-S18
      6  897  911 (S16E05 L=226)   334    165        0 ( 854/ 895- 978/ 941)  23.2 E07-W00 S17-S14
      7 1534  813 (S07W35 L=267)   590    116        0 (1494/ 778-1574/ 874)  80.5 W32-W37 S09-S03
      8  939  829 (S20E01 L=231)   260     67        0 ( 901/ 804- 981/ 847)  45.0 E02-W01 S21-S18
      9  721  730 (S30E13 L=219)   682    181        1 ( 638/ 704- 815/ 768)   4.8 E18-E06 S31-S28
     10 1742  585 (S14W61 L=293)  1796    273        3 (1687/ 530-1784/ 661)  77.5 W56-W66 S17-S08
     11 1657  591 (S16W53 L=284)   878    196        3 (1611/ 535-1695/ 642) -24.0 W46-W56 S19-S13
     12 1444  570 (S23W38 L=270)  1734    348        0 (1287/ 523-1605/ 626)  37.9 W29-W47 S29-S15
     13 1407  350 (S37W47 L=279)   326     64        1 (1379/ 330-1431/ 378)  90.0 W45-W49 S38-S34
  

    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