Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  621 1734 (N27E48 L= 86)   451    123        1 ( 590/1692- 652/1784) -76.0 E50-E46 N23-N31
      2  965 1624 (N27E20 L=115)   990    203        1 ( 875/1559-1050/1678)  -7.6 E27-E12 N24-N29
      3  581 1568 (N15E42 L= 93)   296    103        0 ( 563/1523- 599/1603) -59.0 E44-E39 N12-N17
      4  693 1416 (N08E29 L=106)   438    165        0 ( 617/1356- 761/1447)  45.0 E32-E24 N02-N10
      5  928 1247 (N03E11 L=124)   203     76        1 ( 898/1210- 956/1286) -31.0 E13-E08 N01-N04
      6 1416 1025 (N03W22 L=157)   308     98        0 (1369/ 990-1467/1048)  38.7 W20-W25 N00-N05
      7  229  963 (S27E54 L= 81)   442     89        1 ( 199/ 942- 261/ 980)   8.1 E57-E50 S28-S25
      8  260  631 (S48E53 L= 82)  1203    171        1 ( 214/ 605- 314/ 659)  -5.2 E58-E46 S49-S46
      9 1654  575 (S14W54 L=189)   337     84        0 (1613/ 559-1691/ 592)  33.7 W51-W57 S16-S12
     10 1163  549 (S30W22 L=157)   225     44        0 (1133/ 529-1188/ 568)   0.0 W20-W24 S31-S29
     11  869  207 (S65W26 L=161)   437     69        1 ( 848/ 188- 890/ 231) -29.7 W22-W29 S67-S63
  

    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