Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  949 1524 (N36W11 L=145)   495    128        0 ( 888/1508-1013/1538) -35.5 W07-W14 N33-N38
      2  192 1285 (N42E64 L= 70)   542     83        2 ( 154/1263- 233/1304)  -3.4 E73-E56 N40-N43
      3 1550  894 (S16W27 L=161)   927    330        0 (1425/ 804-1656/1020)  14.0 W17-W37 S18-S12
      4 1183  836 (S09W04 L=138)   771    259        2 (1121/ 717-1235/ 944)  82.4 W02-W05 S16-S01
      5  235  878 (N16E58 L= 76)   596    106        2 ( 199/ 826- 266/ 911) -80.5 E62-E55 N12-N18
      6 1571  791 (S23W27 L=160)   314    121        0 (1537/ 744-1606/ 826)  26.6 W23-W29 S25-S21
      7 1591  697 (S30W26 L=160)   248     53        0 (1563/ 683-1618/ 716) -45.0 W24-W28 S31-S27
      8 1076  634 (S17E08 L=126)   307    122        1 (1030/ 617-1134/ 659) -15.9 E10-E03 S18-S14
      9 1390  568 (S31W09 L=143)   321     45        1 (1366/ 544-1418/ 591) -63.4 W07-W10 S33-S29
     10 1263  490 (S32E01 L=133)   595    124        0 (1218/ 438-1292/ 540) -90.0 E03-W00 S36-S27
     11 1237  398 (S38E06 L=128)   391    115        0 (1179/ 373-1271/ 427)   0.0 E10-E03 S39-S36
  

    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