Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  501 1341 (N38E26 L=163)   703    175        2 ( 442/1273- 571/1400)  18.4 E31-E19 N34-N39
      2  251 1295 (N41E50 L=140)   365     85        1 ( 221/1271- 277/1311)  21.8 E52-E47 N39-N42
      3  871 1230 (N21E03 L=186)   557    163        2 ( 798/1183- 954/1268) -56.3 E07-W00 N16-N25
      4  144 1038 (N27E63 L=127)   562    175        1 ( 108/1021- 189/1051)   0.0 E67-E57 N25-N27
      5  217  962 (N21E56 L=133)  1129    154        2 ( 186/ 921- 264/1030) -66.8 E59-E52 N18-N25
      6  455  834 (N09E39 L=150)   360    106        3 ( 407/ 804- 497/ 879) -50.2 E41-E36 N06-N12
      7  772  772 (S02E21 L=169)   764    172        3 ( 680/ 718- 849/ 811) -36.9 E26-E17 S06-N01
      8  152  750 (N09E73 L=116)  1020    141        2 ( 134/ 703- 172/ 806)  45.0 E78-E69 N06-N12
      9 1661  672 (S34W32 L=221)  1471    183        3 (1542/ 574-1739/ 775)  19.4 W20-W39 S37-S29
     10  429  688 (N01E46 L=144)   304     86        1 ( 406/ 650- 450/ 726)  26.6 E48-E42 N00-N03
     11  944  610 (S15E16 L=173)   484    118        2 ( 913/ 571-1001/ 648)   9.5 E18-E12 S17-S12
     12  684  512 (S15E35 L=155)   576    185        0 ( 655/ 442- 721/ 611)  65.6 E36-E31 S19-S08
     13  448  427 (S15E57 L=132)   668    102        2 ( 401/ 393- 492/ 458) -51.3 E60-E54 S17-S12
  

    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