Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  988 1853 (N58W40 L=336)  2576    348        3 ( 865/1736-1099/1954)   6.4 W17-W80 N53-N61
      2 1529 1721 (N28W65 L=  1)  1620    178        1 (1444/1673-1597/1746) -15.5 W57-W75 N24-N31
      3  542 1602 (N56E21 L=275)   708    130        1 ( 479/1563- 601/1634)   0.0 E28-E13 N53-N57
      4 1451 1420 (N15W37 L=333)   332     65        0 (1427/1383-1468/1452)  78.7 W35-W38 N12-N17
      5  481 1354 (N39E28 L=267)   365    101        0 ( 435/1339- 531/1366)  -6.3 E32-E23 N37-N40
      6 1706  863 (S23W39 L=335)   801    266        2 (1631/ 758-1762/ 974)  32.7 W30-W45 S27-S18
      7 1647  887 (S20W34 L=330)   425    168        0 (1614/ 814-1689/ 952)  32.0 W30-W38 S22-S17
      8 1389  648 (S26W11 L=307)   265     51        0 (1361/ 623-1414/ 673) -68.2 W09-W12 S28-S23
      9 1196  570 (S25E03 L=293)   665    268        0 (1109/ 475-1303/ 624)  14.9 E10-W05 S28-S23
     10 1518  495 (S41W18 L=314)   698    186        1 (1499/ 401-1535/ 579)  74.1 W15-W19 S48-S34
     11  899  486 (S22E23 L=273)   359     97        0 ( 872/ 439- 936/ 517)  38.7 E25-E19 S24-S20
     12 1078  388 (S33E16 L=279)   430    125        0 (1020/ 376-1141/ 404) -12.5 E20-E11 S35-S31
     13  731  351 (S26E41 L=255)   543     95        1 ( 713/ 304- 751/ 402)  36.9 E44-E36 S29-S22
  

    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