Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1591 1601 (N32W49 L=152)  1423    212        3 (1501/1536-1668/1674) -48.0 W43-W54 N26-N37
      2 1403 1554 (N29W30 L=134)  1571    362        2 (1253/1459-1522/1666) -48.8 W21-W37 N22-N38
      3 1020 1623 (N36W04 L=107)   460     80        1 ( 986/1600-1064/1647) -31.0 W01-W06 N34-N38
      4  693 1487 (N29E20 L= 83)   443    149        1 ( 635/1447- 758/1520)  26.6 E24-E15 N26-N30
      5 1012 1397 (N21W01 L=105)   252    117        0 ( 990/1350-1038/1430) -78.7 W00-W02 N17-N22
      6  269 1357 (N23E55 L= 48)   843    100        2 ( 230/1317- 303/1406)  39.8 E59-E51 N20-N26
      7 1080 1343 (N17W05 L=108)   329    132        0 (1016/1323-1136/1371)   0.0 W01-W08 N15-N18
      8 1344 1313 (N14W22 L=125)   195     83        0 (1319/1283-1377/1348) -59.0 W20-W23 N11-N16
      9 1037  871 (S11E00 L=103)   370    127        0 ( 995/ 834-1081/ 929)  45.0 E02-W02 S13-S08
     10 1542  730 (S22W33 L=137)   688    209        0 (1497/ 630-1575/ 833)  80.5 W30-W35 S28-S16
     11 1823  780 (S20W59 L=162)   358     57        2 (1811/ 750-1839/ 814) -76.0 W57-W61 S21-S17
     12 1045  674 (S23E01 L=103)   396    131        0 (1030/ 610-1063/ 736)  82.9 E01-W00 S27-S19
     13 1712  650 (S28W50 L=154)   434     35        2 (1689/ 622-1735/ 678) -45.0 W48-W52 S29-S26
     14 1779  579 (S33W62 L=165)  1691    205        3 (1698/ 493-1849/ 653)  32.7 W54-W69 S38-S28
     15 1540  482 (S39W39 L=142)   833    131        2 (1512/ 425-1563/ 542) -82.9 W35-W43 S43-S34
     16 1656  469 (S41W53 L=156)   534     55        1 (1628/ 445-1690/ 488) -11.3 W48-W58 S42-S39
     17  623  284 (S48E47 L= 56)   670    100        0 ( 598/ 258- 651/ 317)  78.7 E51-E42 S50-S45
     18 1375  246 (S59W34 L=138)   640     91        1 (1344/ 225-1409/ 272) -14.0 W28-W40 S61-S55
  

    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