Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1687 (N51W16 L=290)   562    111        2 ( 900/1648- 965/1728)  16.7 W11-W21 N47-N53
      2 1260 1644 (N37W37 L=311)   347     60        0 (1244/1611-1280/1672)  53.1 W34-W39 N34-N38
      3 1645 1554 (N19W61 L=335)  2240    180        4 (1571/1497-1715/1600) -60.9 W56-W70 N14-N24
      4 1479 1514 (N23W44 L=318)   672    170        1 (1454/1433-1499/1587)  59.0 W40-W48 N17-N27
      5  394 1323 (N39E40 L=234)  1253    172        4 ( 321/1263- 462/1358) -31.0 E49-E33 N33-N42
      6 1625 1321 (N08W46 L=320)   281      0       -1 (1608/1297-1646/1350) -90.0 W44-W47 N06-N09
      7  848 1244 (N24E05 L=269)  1117    425        6 ( 674/1215-1039/1272) -20.0 E17-W05 N18-N27
      8 1628 1096 (S04W39 L=313)   969    145        2 (1590/1028-1654/1171)  49.4 W34-W42 S08-N00
      9 1180 1120 (N09W11 L=285)   409    116        0 (1131/1086-1219/1162)  15.9 W07-W14 N07-N10
     10  675 1092 (N19E20 L=254)   401    120        1 ( 634/1044- 701/1128)  45.0 E22-E17 N16-N20
     11  574 1099 (N21E26 L=248)   257     85        0 ( 527/1079- 610/1116)  -8.1 E30-E23 N20-N22
     12  377 1045 (N22E43 L=231)  1103    180        3 ( 311/ 979- 441/1101) -48.0 E48-E38 N16-N26
     13 1662  803 (S23W36 L=310)  1523    232        0 (1587/ 695-1732/ 894)  26.6 W28-W42 S27-S18
     14  192  803 (N10E68 L=206)   840    118        3 ( 169/ 753- 213/ 869)  66.8 E73-E64 N06-N13
     15 1060  828 (S05E03 L=271)   233     72        1 (1034/ 798-1089/ 864) -68.2 E03-E01 S07-S02
     16 1401  398 (S43W09 L=283)  4204    695       13 (1098/ 204-1550/ 584)  19.6 E23-W23 S51-S32
     17 1046  518 (S23E12 L=262)   289    102        0 (1000/ 486-1089/ 549)  15.9 E15-E08 S24-S22
     18  917  322 (S34E28 L=246)   695    150        0 ( 860/ 284- 977/ 373)  17.1 E34-E21 S36-S31
  

    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