Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1722 1340 (N18W54 L= 34)   723     91        4 (1691/1305-1758/1380) -90.0 W50-W58 N16-N21
      2  212 1250 (N20E64 L=276)   773    110        1 ( 190/1201- 244/1305) -33.7 E68-E58 N17-N23
      3 1606 1226 (N13W41 L= 20)   259     60        0 (1569/1216-1639/1238)  -9.5 W37-W43 N11-N13
      4  408 1072 (N10E41 L=299)   239    103        0 ( 395/1027- 440/1106)  59.0 E42-E38 N06-N11
      5  564 1002 (N05E29 L=310)   698    194        1 ( 539/ 916- 588/1099)  79.7 E31-E27 N00-N11
      6  286  870 (S03E53 L=287)  3500    607        6 ( 209/ 659- 485/1075) -50.2 E61-E38 S16-N09
      7 1611  872 (S09W38 L= 18)   253     86        0 (1574/ 838-1648/ 900)  31.0 W35-W40 S10-S07
      8 1008  865 (S05E02 L=338)   292     92        0 ( 965/ 842-1052/ 898)  21.8 E04-W01 S06-S03
      9 1228  732 (S15W11 L=351)   406    121        2 (1167/ 697-1268/ 760) -33.7 W07-W13 S17-S12
     10  454  581 (S21E44 L=296)   781    174        0 ( 376/ 544- 530/ 624) -24.4 E49-E38 S23-S18
     11  706  428 (S33E28 L=312)   344     32        0 ( 682/ 406- 734/ 448)  14.0 E29-E25 S34-S31
     12  610  299 (S44E44 L=295)  1371    146        1 ( 557/ 240- 652/ 352)  22.9 E62-E36 S50-S38
  

    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