Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  849 1868 (N72W26 L= 72)   674     72        1 ( 824/1835- 870/1900)  11.8 W17-W41 N67-N75
      2 1240 1805 (N50W51 L= 97)   478     62        0 (1217/1773-1267/1831)   8.7 W44-W57 N48-N50
      3 1430 1624 (N32W48 L= 94)   296     78        0 (1419/1583-1443/1661)  29.7 W44-W51 N29-N34
      4  906 1326 (N27E01 L= 45)   675    238        1 ( 854/1235- 962/1419)  42.0 E05-W05 N22-N32
      5  525 1232 (N29E29 L= 17)   587    191        1 ( 493/1139- 545/1316) -80.5 E31-E27 N23-N35
      6 1728 1152 (S03W49 L= 95)  1258    242        2 (1662/1031-1777/1226) -85.2 W44-W51 S10-N02
      7 1603 1085 (S03W37 L= 83)   340     83        0 (1565/1060-1642/1110) -38.7 W34-W39 S05-N00
      8 1851 1071 (S11W60 L=106)   493     82        2 (1829/1027-1868/1107)  38.7 W56-W62 S13-S09
      9  655 1063 (N17E22 L= 24)   280    111        0 ( 625/1022- 694/1101)  33.7 E24-E18 N14-N18
     10  410 1008 (N18E41 L=  5)   296     74        1 ( 376/ 982- 437/1043) -51.3 E43-E38 N15-N20
     11 1407  775 (S16W18 L= 64)   340    146        1 (1353/ 749-1479/ 817)   0.0 W14-W23 S16-S14
     12 1686  657 (S32W39 L= 85)   430    102        1 (1646/ 640-1737/ 678)  -6.3 W34-W44 S33-S30
     13 1316  585 (S26W08 L= 54)   429     98        1 (1292/ 548-1336/ 648) -90.0 W06-W09 S28-S21
     14  513  440 (S19E50 L=356)   372     57        1 ( 496/ 410- 529/ 469)  38.7 E53-E48 S21-S17
     15  419  423 (S19E61 L=345)   501    102        2 ( 376/ 412- 468/ 437) -18.4 E66-E56 S20-S17
     16  698  374 (S27E39 L=  7)  1343    245        2 ( 601/ 328- 817/ 430)   8.1 E49-E28 S29-S24
  

    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