Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  535 1683 (N25E54 L=107)  1460    304        0 ( 496/1618- 578/1758) -39.8 E62-E49 N19-N29
      2  368 1566 (N14E62 L= 99)  1938    385       10 ( 330/1435- 411/1663) -63.4 E69-E56 N04-N20
      3  417 1591 (N16E59 L=102)   534    111        2 ( 396/1562- 433/1620) -76.0 E60-E56 N14-N18
      4  551 1417 (N08E39 L=121)   448    102        1 ( 528/1358- 574/1451) -63.4 E41-E37 N04-N10
      5  906 1173 (N02E10 L=150)   305     74        1 ( 871/1129- 941/1214)  63.4 E11-E08 N00-N05
      6 1376  965 (N02W22 L=182)   411    139        1 (1313/ 899-1419/1020)  69.4 W19-W23 S03-N05
      7 1365  601 (S19W32 L=193)   339    100        1 (1313/ 581-1407/ 636)  31.0 W29-W34 S20-S16
      8 1592  571 (S14W50 L=210)   621    117        0 (1542/ 526-1637/ 624)  76.0 W47-W51 S18-S10
      9 1681  571 (S12W57 L=218)   544     97        1 (1639/ 542-1718/ 595)  33.7 W54-W60 S13-S09
     10 1037  430 (S38W19 L=180)  1142    324        1 ( 911/ 405-1213/ 464)  18.4 W09-W30 S42-S32
     11 1435  313 (S34W53 L=214)   520    163        0 (1399/ 295-1499/ 342)  18.4 W48-W57 S35-S30
  

    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