Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1067 1853 (N55W47 L=272)   902    103        0 (1032/1813-1108/1884)   6.3 W38-W56 N53-N56
      2 1698 1566 (N15W67 L=292)  2069    191        4 (1618/1491-1758/1626) -59.0 W62-W72 N09-N19
      3  886 1453 (N34W04 L=229)  1527    444        3 ( 699/1409-1104/1486) -30.3 E08-W16 N25-N39
      4 1618 1399 (N09W48 L=273)   560    112        1 (1579/1354-1647/1452)  24.0 W44-W53 N07-N11
      5 1094 1249 (N17W10 L=235)   551    248        0 ( 970/1221-1195/1289) -12.1 W02-W16 N14-N19
      6 1284 1210 (N09W20 L=244)   629    158        1 (1210/1157-1339/1262) -66.0 W16-W22 N04-N13
      7  729 1167 (N22E14 L=211)   733    219        2 ( 634/1130- 826/1196)  -4.4 E20-E07 N19-N23
      8  255 1137 (N31E53 L=172)   487     73        0 ( 239/1097- 269/1178) -68.2 E55-E51 N28-N33
      9 1812 1079 (S14W53 L=278)   369     64        1 (1789/1050-1840/1100)  11.3 W51-W57 S15-S12
     10  735 1040 (N14E16 L=209)   594    270        0 ( 648/ 986- 830/1072) -15.3 E22-E11 N10-N16
     11  232  892 (N17E59 L=166)   961    211        0 ( 156/ 837- 318/ 933) -24.0 E68-E50 N12-N20
     12 1549  804 (S21W25 L=250)   352    102        0 (1516/ 761-1583/ 845)  26.6 W22-W28 S22-S19
     13 1737  582 (S44W42 L=267)  2803    317        4 (1628/ 496-1825/ 716)  22.2 W28-W52 S47-S37
     14  739  454 (S19E34 L=190)   378     86        1 ( 722/ 409- 752/ 498)  50.2 E36-E31 S22-S16
     15 1458  387 (S48W11 L=235)   806    125        1 (1398/ 357-1518/ 416)  -4.8 W04-W16 S50-S45
     16 1161  212 (S51E23 L=202)  1223    169        3 (1079/ 182-1243/ 240) -10.0 E31-E14 S54-S48
  

    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