Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  368 1603 (N12E64 L=184)  2570    385        7 ( 259/1518- 459/1695)  49.4 E73-E58 N05-N19
      2 1224 1332 (N15W03 L=251)   725    203        0 (1162/1250-1281/1427) -35.0 E02-W08 N12-N19
      3  568 1249 (S06E32 L=215)   233     82        0 ( 525/1240- 608/1258)  18.4 E35-E29 S06-S04
      4  432 1063 (S19E37 L=210)   289     94        1 ( 394/1027- 454/1089)   0.0 E40-E35 S20-S17
      5  607 1068 (S15E25 L=223)   228     72        1 ( 559/1055- 642/1086)  31.0 E28-E23 S16-S13
      6  383  879 (S31E38 L=210)   630    121        1 ( 337/ 857- 441/ 909)   0.0 E42-E33 S31-S29
      7 1019  819 (S18W05 L=253)  1040    312        1 ( 866/ 780-1183/ 848)  18.4 E03-W15 S22-S15
      8  706  817 (S27E13 L=235)   485    112        2 ( 646/ 801- 764/ 837)  14.0 E17-E09 S28-S24
      9 1297  632 (S21W28 L=275)  1393    414        1 (1129/ 533-1466/ 703)  11.3 W17-W42 S25-S16
     10 1602  541 (S16W52 L=300)   395     43        1 (1578/ 516-1626/ 563)  71.6 W50-W54 S17-S14
     11 1603  437 (S21W58 L=306)   493    143        0 (1585/ 400-1619/ 474) -14.0 W54-W62 S22-S19
     12 1336  323 (S36W47 L=295)   366     42        1 (1314/ 303-1362/ 343)  90.0 W44-W49 S37-S34
  

    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