Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1689 1629 (N19W69 L=173)  1669    140        2 (1666/1559-1721/1683)  14.7 W63-W83 N15-N21
      2 1302 1666 (N33W39 L=143)   455     88        1 (1253/1645-1350/1680) -33.7 W36-W42 N30-N35
      3 1324 1570 (N27W35 L=139)   778    209        1 (1248/1500-1407/1606)   7.6 W27-W42 N24-N28
      4 1393 1418 (N16W32 L=136)   276     88        0 (1361/1379-1411/1454) -90.0 W30-W33 N14-N19
      5 1125 1351 (N19W14 L=118)   247     74        0 (1093/1314-1151/1386) -80.5 W13-W15 N16-N22
      6  850 1280 (N22E04 L=100)   850    220        2 ( 795/1197- 933/1354) -71.6 E06-W00 N15-N27
      7 1738 1160 (S07W48 L=152)  2293    414        1 (1683/ 967-1779/1322)  63.4 W40-W52 S16-N02
      8  597 1215 (N25E21 L= 83)  2145    543        4 ( 365/1080- 794/1293) -34.2 E39-E10 N12-N32
      9  480 1081 (N19E32 L= 72)   227    110        0 ( 427/1065- 535/1097) -29.7 E35-E28 N17-N21
     10 1842  861 (S27W53 L=157)  2676    243        4 (1727/ 800-1934/ 931)   4.8 W40-W68 S29-S24
     11 1545  844 (S19W27 L=130)  1581    330        3 (1418/ 744-1670/ 917)   8.5 W16-W36 S21-S16
     12  363  602 (S05E54 L= 50)  7051      0       21 ( 251/ 446- 461/ 831)  55.7 E65-E43 S15-N07
     13  750  719 (S07E23 L= 81)   355    129        0 ( 718/ 672- 770/ 769)  51.3 E26-E20 S08-S03
     14 1604  710 (S29W29 L=133)   307     72        0 (1568/ 687-1629/ 741)   9.5 W25-W31 S30-S27
     15  515  352 (S23E56 L= 47)   494     65        0 ( 481/ 330- 540/ 379) -76.0 E58-E53 S25-S20
     16  630  337 (S27E48 L= 56)   708    128        2 ( 572/ 311- 687/ 376)   0.0 E54-E41 S28-S24
     17 1410  233 (S62W05 L=109)   464     45        0 (1386/ 215-1440/ 254) -48.8 W02-W09 S66-S58
  

    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