Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  648 1638 (N30E34 L=193)   390     53        0 ( 616/1620- 684/1652) -11.3 E36-E31 N28-N30
      2  404 1379 (N11E45 L=182)   324    113        0 ( 392/1319- 417/1426) -71.6 E46-E44 N07-N13
      3  480 1247 (N03E37 L=190)   468     70        0 ( 456/1200- 506/1294)  78.7 E37-E35 N00-N06
      4  713 1143 (S01E20 L=207)   185     21        0 ( 697/1121- 730/1165)  63.4 E20-E18 S02-N00
      5  378  961 (S14E42 L=185)   250     61        1 ( 354/ 934- 394/ 996)  76.0 E43-E40 S16-S12
      6  843  918 (S13E10 L=217)   279     64        0 ( 806/ 899- 881/ 938)  21.8 E12-E07 S14-S12
      7 1408  796 (S15W27 L=253)  1359    274        2 (1278/ 716-1531/ 882)  40.6 W19-W33 S20-S07
      8  658  773 (S24E21 L=206)   263     95        0 ( 620/ 751- 711/ 784)  18.4 E23-E17 S25-S22
      9  835  697 (S27E09 L=218)   544    173        0 ( 761/ 642- 903/ 748)  41.6 E13-E04 S31-S23
     10 1692  682 (S17W50 L=277)   281     53        0 (1671/ 648-1713/ 706)  76.0 W48-W51 S19-S15
     11 1584  637 (S21W43 L=270)   359     87        0 (1537/ 595-1614/ 677)  63.4 W40-W44 S24-S18
     12  660  344 (S55E25 L=201)   474     81        1 ( 634/ 318- 697/ 370)  20.6 E29-E21 S57-S52
  

    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