Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  488 1531 (N22E42 L=212)  1695    627        3 ( 344/1448- 666/1655) -20.9 E61-E27 N17-N30
      2  353 1147 (S03E45 L=209)   512    186        2 ( 293/1110- 429/1173)  -9.5 E50-E38 S04-S01
      3  421 1033 (S09E39 L=216)   224     72        0 ( 392/ 994- 446/1068)  59.0 E40-E37 S11-S06
      4  185  984 (S13E60 L=194)   439    110        3 ( 155/ 960- 217/1011)  -7.1 E64-E56 S14-S11
      5  395  950 (S14E41 L=214)   285    120        0 ( 346/ 929- 435/ 979)  15.9 E44-E37 S15-S12
      6  967  875 (S15E02 L=252)  1098    254        2 ( 857/ 791-1099/ 965)  40.2 E08-W05 S20-S08
      7  740  941 (S13E17 L=238)   227     73        1 ( 699/ 918- 777/ 958)   0.0 E19-E14 S13-S11
      8  514  921 (S15E32 L=223)   250     62        1 ( 479/ 899- 545/ 938) -11.3 E34-E29 S16-S14
      9  196  903 (S18E59 L=196)   739    151        4 ( 166/ 861- 226/ 936) -51.3 E62-E55 S20-S15
     10  650  862 (S18E22 L=233)   232     81        1 ( 632/ 834- 676/ 908)  76.0 E23-E20 S19-S15
     11 1230  749 (S20W16 L=270)   406    132        0 (1197/ 695-1279/ 816)  69.4 W14-W17 S23-S15
     12  745  788 (S22E15 L=239)   370     94        1 ( 695/ 756- 795/ 805)   0.0 E18-E12 S24-S20
     13  296  763 (S26E50 L=204)   411    115        0 ( 279/ 715- 327/ 804)  50.2 E51-E46 S29-S23
     14  439  723 (S28E38 L=217)   559    186        0 ( 401/ 654- 494/ 774)  45.0 E41-E33 S32-S24
     15 1595  595 (S24W45 L=299)  1139    158        1 (1522/ 537-1659/ 649)  48.8 W40-W48 S28-S20
     16 1481  530 (S30W38 L=293)   321    130        1 (1463/ 471-1506/ 587) -74.1 W36-W40 S33-S26
     17 1661  397 (S35W60 L=314)   549    137        3 (1625/ 369-1689/ 416)  21.8 W56-W62 S37-S33
  

    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