Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1189 1693 (N48W24 L= 48)   576     88        1 (1158/1645-1214/1727) -90.0 W21-W26 N43-N51
      2 1588 1421 (N23W46 L= 69)   861    170        1 (1526/1354-1643/1479) -60.9 W42-W48 N18-N27
      3  257 1414 (N34E66 L=318)  1106    105        2 ( 227/1381- 301/1451) -10.0 E74-E57 N31-N36
      4 1431 1342 (N20W31 L= 54)   295    111        1 (1391/1304-1477/1372)  23.2 W27-W34 N18-N22
      5 1127 1058 (N06W07 L= 30)   229     53        1 (1095/1044-1155/1077)  14.0 W04-W08 N05-N07
      6 1350 1022 (N02W20 L= 44)   372    131        1 (1311/ 973-1389/1075)  38.7 W17-W23 N00-N04
      7 1858  947 (S10W62 L= 86)  2069    214        3 (1821/ 854-1888/1028)  68.2 W57-W66 S15-S04
      8 1917  939 (S11W72 L= 96)  1069     83        2 (1895/ 899-1934/ 976)  78.7 W67-W76 S13-S08
      9  309  839 (S02E52 L=332)   455     95        0 ( 278/ 800- 334/ 882) -45.0 E54-E49 S04-N00
     10  245  717 (S10E62 L=322)  2340    336        3 ( 194/ 584- 319/ 812) -65.0 E66-E54 S18-S03
     11 1496  741 (S17W28 L= 52)   230     72        0 (1479/ 711-1520/ 779)  53.1 W27-W30 S19-S14
     12  410  698 (S11E45 L=338)   336     53        2 ( 382/ 673- 433/ 728)  18.4 E48-E42 S12-S09
     13  861  615 (S19E14 L=  9)   530    181        0 ( 795/ 571- 896/ 684)  40.6 E19-E11 S21-S14
     14  945  441 (S32E12 L= 12)  3952    892        2 ( 722/ 221-1327/ 675)  26.1 E39-W16 S50-S19
     15 1601  608 (S28W38 L= 62)   331     79        0 (1562/ 590-1638/ 628)   0.0 W34-W41 S28-S26
     16  590  478 (S27E37 L=346)  1975    522        3 ( 488/ 328- 678/ 615) -76.6 E46-E31 S38-S16
     17 1689  531 (S35W50 L= 74)  1164    121        2 (1620/ 504-1747/ 558)  -7.1 W42-W59 S36-S32
     18 1043  395 (S37E05 L= 18)   375    111        1 ( 988/ 384-1101/ 412)  -6.3 E09-W00 S37-S35
  

    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