Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1704 (N35E23 L=196)   732    100        1 ( 823/1664- 912/1731)  51.3 E25-E19 N31-N37
      2  969 1598 (N29E12 L=207)   249     73        0 ( 941/1562-1001/1617)  68.2 E13-E10 N25-N30
      3 1173 1399 (N18W04 L=223)   251     23        0 (1153/1379-1194/1425)  90.0 W02-W05 N17-N20
      4 1621 1214 (N14W35 L=254)   287     68        1 (1594/1195-1651/1244)  14.0 W33-W37 N13-N15
      5  784 1145 (S02E16 L=203)   395    156        1 ( 760/1073- 814/1230) -60.3 E17-E13 S05-N03
      6 1917 1140 (N16W64 L=283)  2079    219        3 (1866/1083-1952/1191) -35.5 W57-W69 N13-N19
      7  584  974 (S14E26 L=192)  1065    342        1 ( 461/ 865- 716/1089)  45.0 E34-E19 S21-S05
      8  380  985 (S15E41 L=177)   771    249        1 ( 284/ 936- 471/1032)  21.8 E49-E34 S18-S11
      9 1219  811 (S16W15 L=233)   497    167        2 (1155/ 752-1262/ 857) -26.6 W10-W18 S18-S13
     10 1269  731 (S20W19 L=238)   195     75        0 (1221/ 715-1302/ 748)  21.8 W16-W21 S21-S18
     11  518  646 (S35E31 L=188)   302     50        1 ( 489/ 626- 538/ 666)   0.0 E33-E29 S36-S33
     12  631  453 (S47E23 L=196)   261     54        0 ( 602/ 440- 663/ 466)   0.0 E25-E20 S48-S46
     13  862  213 (S65W02 L=221)   540     63        1 ( 835/ 195- 894/ 240)  -7.1 E01-W07 S67-S62
  

    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